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We propose coalescent mechanism of economic grow because of redistribution of external resources. 
It leads to Zipf distribution of firms over their sizes, turning to stretched exponent because of size- 
dependent effects, and predicts exponential distribution of income between individuals. 

We present new approach to describe fluctuations on the market, based on separation of hot (short- 
time) and cold (long-time) degrees of freedoms, which predicts tent-like distribution of fluctuations 
with stable tail exponent (x = 3 {(X = 2 for news). The theory predicts observable asymmetry 
of the distribution, and its size dependence. For financial markets the theory explains first time 
"market mill" patterns, conditional distribution, "D-smile", z-shaped response, "conditional double 
dynamics" , the skewness and so on. 

We propose a set of Langeven equations for the market, and derive equations for multifractal 
random walk model. We find logarithmic dependence of price shift on the volume, and volatility 
patterns after jumps. We calculate correlation functions and Hurst exponents at different time 
scales. We show, that price experiences fractional Brownian motion with chaotically switching of 
sub- and super-diffusion, and calculate corresponding probabilities, response functions, and risks. 
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I. INTRODUCTION 

First question behind any research is why do we need 
it? There are no unique approach in econophysics, and 
the number of different approaches grows exponentially 
with time. How can we decide, which of them is "cor- 
rect" , if by construction, any one well describes empirical 
facts? 

The answer is simple: in no way. All of them are 
equivalent at regions of their applicability. But these re- 
gions are very different, and only several theories able to 
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describe a large variety of empirical facts. Extrapolat- 
ing, one can assume that only one theory can predict all 
important phenomena: Market mill patterns, multifrac- 
tality of fluctuations, volatility patterns, different Hurst 
exponents above and below a time t x , and many other 
facts. Some of them can be described in different ways, 
but not the hookup of all facts. 

What criteria should satisfy such theory? At first sight, 
it is mathematical rigor. The most striking example is 
the Flory approach in polymer physics, which is abso- 
lutely "wrong" mathematically, but extremely well de- 
scribing all known situations. All multiple attempts to 
(im)prove it were failed. We conclude, that the rigor 
of the theory is usually "inversely proportional" to intu- 
ition. 

Well, what kind of the theory should not be? If any 
new fact or their series need to introduce additional terms 
or ideas into the theory, the later can be considered as 
a collection of facts, arbitrary ordered according to the 
test of the author. We think, the real theory must predict 
in future yet unknown facts (at present this criterium is 
equivalent to extremely wide region of its applicability), 
to be as rigor as possible, and use minimum of initial 
assumptions. 

We do not know other criteria of the "validity" of the 
theory, and this is the reason why the theory must de- 
scribe all known trustable facts. Present paper can be 
considered as an attempt to follow this criteria. Only 
one main idea lays in the basis of our theory of market 
fluctuations: we assume, that they can be described as 
random walk motion at all time scales. In the case of fi- 
nancial market, it is random trading at all time horizons 
from seconds to tenths years. 

Our theory can be considered as an attempt to make 
a step from numerous descriptive approaches toward a 
physical Langeven formulation of the "econophysical" 
problem. This is why we emphasize analogies with other 
branches of physics, which may confuse econo-physicists 
otherwise. Although we show, that multi-time random 
trading allows to explain most of market dynamics, it 
may be extended later in many directions. 

As a strategy line, for each problem we try to construct 
a simplified model of such multi-time random motion, 
capturing the most of physics. As the result, we left 
with several parts of the whole puzzle, strongly inter- 
correlated with each other. It is the reason of unusual 
length of this paper, which can not be cut into several 
independent small parts. 



II. FIRMS, CITIES AND INCOME 
DISTRIBUTIONS 

A. Is there thermodynamics of the market? 

Econophysics studies physical problems in economics, 
and most of its results were obtained from analogy with 
thermodynamics. One of classical problems of econo- 



physics, the firm grow, is usually described by the model 
of stochastic firm growing 1 . In order to explain empir- 
ically observed Zipf distribution of firm sizes 2 it is pro- 
posed to introduce the lower reflecting boundary in the 
space of firm sizes, which stabilizes the distribution to a 
power law 3 . Unfortunately, this explanation is inconsis- 
tent for firms of one or several employers, well described 
by the same empirical Zipf distribution. 

Different models of internal structure of firms were pro- 
posed for the stochastic mechanism of firm growing. Hi- 
erarchical tree- like model of firm was studied in Refs. 4 ' 5 . 
A model of equiprobable distribution of all partitions of a 
firm was introduced in Rcf. 6 . Both models neglect the ef- 
fect of competition between different firms. The random 
exchange of resources between firms was taken into con- 
sideration in "saving" models 7 . In Refs. 8-11 the process 
of stochastic firm grow and loss was considered by anal- 
ogy with scattering processes in liquids and gases. The 
distribution of firms over their sizes in different countries 
was studied in Rcf. 12 . 

The theory of firms is usually called microeconomics, 
and from economical point of view it is hard to con- 
sider the stochasticity as the moving force of economic 
grow. While in thermodynamics the stochasticity origi- 
nates from interaction with a huge "thermostat", there 
are no such thermostat for the market, which subsists 
only because of activity of its direct participants. 

This puzzle forces us to develop a "mean field" theory 
of firm growing, neglecting any fluctuation processes. We 
show, that the moving force of evolution on the market 
are not thermal- like excitations, but the supply of exter- 
nal resources, which are (re-)distributed between differ- 
ent firms. Exhaustion of the resource kills this (part of 
the) market, while appearance of a new resource gives 
rise to a new market. The process of firm growing and 
mergence is similar to coalescence of droplets of a new 
phase, when stochasticity plays only minor role. 

In section II B we show that the coalescence theory 
predicts Pareto power low for the distribution of firm 
sizes. We propose self-similar tree-like model of firms 
in section II B 2. This model is solved in Appendix B, 
and we show, that it explains empirically observable time 
dependence of the Pareto exponent for the world income. 

The formal resemblance of observable exponential dis- 
tribution of the income between individuals to Boltz- 
mann statistics was used in Ref. 13 to justify the appli- 
cability of methods of equilibrium thermodynamics. But 
how can all sectors of country economics and services 
always be in thermal equilibrium? In section II B 3 we 
propose an alternative explanation, based on unified tax 
policy in the whole country: the coalescent approach pre- 
dicts, as a by-product, the exponential income distribu- 
tion, even without invention of thermal equilibrium. This 
distribution is valid for the majority of the population, 
and statistical fluctuations are only responsible for power 
tails of its upper part (1-3%). 

Countries with different financial policy have different 
"effective temperature" of the distribution, which can 
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be equilibrated only after unification of their financial 
policies, even without establishment of a "heat death" - 
global thermal equilibrium. Although one may consider 
the perpetual trade deficit of US as consequence of the 
fundamental second law of thermodynamics 13 , it would 
be more natural to explain it by financial policy, directed 
on attraction of resources to the country. 

Econophysics is not only one field, deceptively resem- 
bling thermodynamics, we have to mention also a sand, 
turbulence and other macroscopic systems, which form 
complex dissipative structures in the response on some 
external forces. Although such "open systems" can not 
be characterized by thermodynamic potentials, the pro- 
cess of dissipation is accompanied by the rise of infor- 
mation entropy. We calculate the entropy of the market 
and show, that it can only increase with time, since the 
market irreversibly absorbs external information (there 
is deep analogy with physics of decoherence, discussed in 
Conclusion) . 

"Thermodynamic-type" models predict asymptotically 
Gaussian distribution of firm grow rates, while the real 
distribution has tent-like shape. In order to reproduce 
it, in Rcf. 14 an artificial potential was introduced in dif- 
fusion equation, restoring the firm size to a certain ref- 
erence value, at which the grow rate abruptly changes 
its sign. In this paper we elaborate a new approach to 
study dynamics of temporal dissipative structures on the 
market, which do not use these artificial assumptions. 

In section II C we introduce new general approach to 
study market fluctuations. Main ideas of this approach 
will be first formulated for the problem of firm grow. The 
market is the system with multiple (quasi-) equilibrium 
states, characterized by extremely wide spectrum of re- 
laxation times. By analogy with glasses, for given obser- 
vation (coarse graining) time interval r we can divide all 
degrees of freedom of the market into "hot" and "cold" 
ones, depending on their relaxation times. Hot degrees of 
freedom are in equilibrium, and they generate high fre- 
quency fluctuations because of uncertainty on the mar- 
ket, while cold degrees of freedom arc not equilibrated, 
and evolve on times large with respect to r. As in the case 
of spin-glasses, high degeneracy of quasi-equilibriums in 
the market is reflected in the presence of a gauge invari- 
ance. Any averages should be defined in two stages: first, 
the annealed averaging over hot degrees of freedom, and 
then quenched averaging over cold degrees of freedom. 

We demonstrate, that our theory reproduces empiri- 
cally observable (in general, asymmetric) tent-like distri- 
bution of firms over their grow rates. In section II C 4 
we show, that this distribution has fat tail with stable 
exponent /x, equals to the number of essential degrees 
of freedom of the noise (/z = 3 for Markovian statistics 
of hot degrees of freedom, and \x — 2 for uncorrelated 
noise) . 



B. Mean field theory 

Dynamics of firm growing is similar to kinetics of grow- 
ing of droplets of a new phase. Large firms can absorb 
smaller ones, and they can grow or leave the business, by 
analogy with resorption and growing of droplets in the 
supersaturated solution. Below we use this analogy to 
construct a new theory, not relying on stochastic mecha- 
nisms of firm growth. Entropic and microeconomic inter- 
pretations of our theory are discussed in Appendixes A 
and C. 



1. Zipf distribution 

For definiteness sake we define the firm size as the 
number G of its employees. In general, it could be any 
resource, shared between different firms on the market. 
According to economic approach (analog of the mean 
field approach in physics) firms can hire or loose the staff 
only through the "reservoir" of unemployments of value 
U (t) at time t. Diffusion processes lead to finite value 
> of the "natural unemployment" . "Actual unem- 
ployment" U is the sum of U* 15 and the "market unem- 
ployment", A(t): 

U (t) = U* + A (t) . 

The equation of the resource balance can be written in 
the form 



Q(t) = U(t)+ / Gf(G,t)dG 



(1) 



where Q (t) is the supply of external resources. The prob- 
ability distribution function (PDF) / (G, t) of firm sizes 
is determined by the continuity equation, 



df(G,t) 

at 



d_ 

dG 



(2) 



where dG/dt is the rate of ordered motion in the space of 
firm sizes. Diffusion contribution in Eq. (2) is negligible 
in coalescent regime. According to the famous Gibrat's 
observation 1 the relative grow rate of the firm, 



I_dG 
G~dt 



(3) 



do not depend on its size, G. In the case of full employ- 
ment, A = 0, the average numbers of people getting a job 
and leaving it are the same, and there are no source for 
firm grow, tq = 0. At small A we can hold only linear 
term in the series expansion of the grow rate tq — qA in 
powers of A with constant q. 

To solve the set of equations (I) - (3) we substitute 
Eq. (3) with tq — qA into Eq. (2), and find its general 
solution 



f{G,t) = ± X 



ln£- - q £ A (if)dt' 
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where Go is the firm size at initial time t = to and \ 
is arbitrary function. Substituting this solution into the 
balance equation (1) and introducing new variable of in- 
tegration u = In (G/Go), we find 

rt 



•(t) = U (t) + Go / e u X 



u-q A(t')dt' 



du. (4) 



Consider the case of power growing of external re- 
sources, 



Q (t) = Q t r 



(5) 



For general time dependence Q (t) its logarithmic rate m 
is determined by expression 



dlnQ(i) 
dint 



(6) 



In the case of small unemployment value, U <^ Q, gen- 
eral solution of Eq. (4) takes exponential form, x i u ) — 
Xo e ~ KU - Substituting this expression into Eq. (4) and 
taking into account that the distribution / (G, t) can not 
depend on initial firm size, Go, we find k — 1 and 



where G m \ n and G,, 



G n 



G,i 



■ exp 



q I A(t')dt' 
o 



are maximal and minimal firm 



sizes on the market. The solution of this equation has 
the form 

A (t) = m/ (qt) , x = Qo/ In (G max /G min ) . (7) 

First of Eqs. (7) predicts, that the economic grow, see 
Eq. (5), leads to less actual unemployment, A, in qual- 
itative agreement with the famous macroeconomic "Fil- 
lips curve". Close quantitative relation between the co- 
alescent theory and the Fillips low is established in Ap- 
pendix C. 

We conclude, that for any monotonically increasing 
function Q (t) the distribution of firms over their sizes 
G has Zipf form: 



f(G,t) = 



'(*) 



1 



In (G max / G m i n ) G 2 



(8) 



This dependence was really observed for extremely wide 
range of firm sizes, see Fig. 1, where empirically observ- 
able distribution 

J^-f(G,t)dG 1_ 
" G 



F (G, t) 



J°—f(G,t)dG 



(9) 



is plotted. The Zipf distribution 2 (8) is valid for the 
entire range of US firms 16 (from G m i n = 1 to G max = 10 6 ) 
with Pareto exponent very close to unity. 

The same mechanism may be responsible for power 
distribution function of cities over their population, the 
amount of assets under management of mutual funds 17 , 
banks 18 and so on. In the analysis of city population in 
different countries, the exact form of Zipf 's law (9) was 
confirmed in 20 out of 73 countries 19 . Deviations from 
this low will be studied in next section. 
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FIG. 1: Size distribution of U.S. business firms in 1997 (Cen- 
sus data) 16 . Straight line corresponds to power law distribu- 
tion F (G) ~ G~~' with exponent 7 = 1.059. 



2. Stretched exponent 

The Pareto exponent (9) can deviate from 1 because of 
ineffective management, strong influence of industry ef- 
fects on small firms and so on. With increasing size, these 
effects gradually trail off, while remaining international, 
national and regional shocks equally affect all firms. As- 
suming self-similiarity of firm structure, the variation of 
the firm size can be described by Master equation 



= G- 1 dG/dt = r G - P G 



-0 



(10) 



with constant p and (3. 

To derive Eq. (10), consider the firm as the self-similar 
tree 14 of n generations, each of Go ^ 1 branches. The 
size Go of each subdivision is described by the same type 
of equation (10), 



r = G 1 dG /dt = r Q 



(11) 



Substituting the estimation G ~ Gq for the size of the 
whole tree in Eq. (10) and comparing with Eq. (11), 
we find the relation between coefficients of Master equa- 
tions (10) and (11): 

j3 = (3 /n, r G = r n, p = p Q n. 

In Appendix C we show that while the Gibrat grow 
rate r G is fixed by economic factors, the coefficient p of 
job destruction can experience strong random fluctua- 
tions Ap. Neglecting fluctuations of r G in Eq. (10) we 
find that fluctuations in size are inversely correlated to 
the size with an exponent /3: 



Ar = -ApG 



-13 



(12) 



In order to estimate the exponent (3 , consider a hy- 
pothetical structureless firm with n = 1 of the size 
G = Go ^> 1. Fluctuations of its size are characterized 
by Gaussian exponent j3 — /3 Q — 1/2. The exponent (3 of 
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real firms takes small values (3 = 0.15— 0.21 5 , correspond- 
ing to the number of tree generations n = 1/ (2(3) =3 — 4. 
Using Eq. (12) we find the dependence of the standard 
deviation of grow rate Ar on the firm size, 



(Ar 2 ) 1/2 =aG-P 



(13) 



1/2 

where a = (Ap 2 ^ does not depend on firm size G. 
This relation is in excellent agreement with empirical 
data 14 - 20 . 

The condition r = (10) determines the critical firm 
size 



G c = ( P /r G ) 



1//3 



(14) 



Small firms with G < G c collapse with time and may 
leave from the business (or reach a certain fluctuation 
size), while large firms with G > G c grow. In Appendix A 
we find the entropy S (G) of the firm of size G, and show 
that G — G c corresponds to its minimum, and also to 
the minimum point of a "U-shaped" average cost curve 
in the conventional economic theory (Appendix C). We 
also derive maximum entropy principle for the market 
(Appendix A), which is known as the most foundational 
concepts of Gibbs systems. 

In Appendix B we show that the solution of 
Eqs. (1), (2) with the rate (10) has stretched exponent 
form: 



F (G) = exp [- (1//3 - m) (G/G c f 



(15) 



Taking the limit (3 — > we reproduce Eq. (9). It is 
shown that stretched exponent is the best fitting ap- 
proximation for many observable distributions (size of 
cities, population of different countries, popularity of ex- 
ecutors, lifetime of different species, strength of earth- 
quakes, indices of quoting, number of coauthors, rela- 
tive rates of protein synthesis and many others 21-23 ), 
which are determined by the competition of units for 
common resources. At small but finite (3 <C 1 expanding 



(G/G c f ~ 1 + /31n (G/G c ) in Eq. (15) we find 



F (G) ~ GT 



7 = 1— (3m. 



(16) 



We conclude, that the exponent 7 of Pareto distribu- 
tion is, in general, not universal and depends on current 
rate m (t) of external supply, Eq. (6) . This conclusion 
can be verified by empirical observations: typically, the 
value of this exponent is in the interval 0.7 < 7 < 1. 
For example, the size distribution of Danish production 
companies with ten or more employees follows a rank-size 
distribution with exponent 7 = 0.741 . 

To confirm the dependence of the exponent 7 on the 
supply rate m (t) , consider the distribution of world in- 
come across different countries. We assume, that coun- 
tries could be described by the same Master equation (10) 
as large firms. Exponential growing of consumable re- 
sources leads to linear time dependence of m ~ t, see 
Eq. (6). As the result, the exponent 7 linearly decreases 
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FIG. 2: Temporal path of the exponent 7 (continuous line), 
and its approximation by linear dependence (dotted line) 25 . 



with time, in good agreement with empirical observa- 
tions, see Fig. 2. Assuming, that Q (t) doubles every 12 
years, we estimate (3 ~ 0.1, corresponding to a reason- 
able number n ~ 5 of hierarchical management ranks in 
the "typical" country. 



3. Income distribution 

In order to find the distribution of income between 
individuals we first introduce the most important eco- 
nomic terms. The total income per state, Q, is shared 
between all individuals {G} and the state expenses, U (t), 
according to the balance equation (1). There are some 
minimal expenses of the state, [/*, and the inequal- 
ity A = U — U* > is usually regulated indirectly, 
through taxes, which determine the relative income rate, 
Tg = <?A, of individuals. Therefore, the income G can be 
described by a generalization of the Master equation (3), 



dG 
~dt 



r G G - p. 



(17) 



The last term describes the rate of losses (living- wage) , 
the same for all individuals (linear in G losses renormalizc 
rc)- Since Eq. (17) has the form of Eq. (10) with (3 — 
1, from Eq. (15) we get exponential distribution of the 



/(G,i)~e- G / T , T = p/[r G (l-m)] 



(18) 



According to Eq. (B3) of Appendix C the average in- 
come (the "temperature 26 " ) T linearly grows with time, 
in good agreement with empirical observations 26 , and 
also rises with the supply rate m (6). It is small for coun- 
tries with low living wage p, producing high inequality in 
incomes. 

Analysis of empirical data shows 13 , that for approxi- 
mately 95% of the total population, the distribution is 
exponential, while the income of the top 5% individuals 
is described by a power-law (16) with time dependent 
Pareto index 7. This tail is because of speculation in 
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stocks, when the income is proportional to the volume 
of sale/buy G <~ V. The distribution of large volumes is 
power tailed, P (V) <~ V" 1 . The exponent 7 is not uni- 
versal, it depends on individual stocks with typical value 
7 ~ 3/2, in good agreement with observable values 26 
7 = 1.4 — 1.8 (changing of the most profitable stocks 
leads to variations in 7). 

In general, the income may come from different 
sources. In the case of n independent sources convolution 
of n exponential distributions gives the Gamma distribu- 
tion P n (G) <~ G"e _G / T , which better describes Russian 
Rosstat data of salary distribution. 



C. Fluctuation theory 

1. Cold and hot degrees of freedom 

Our approach to the description of fluctuations on the 
market is related to the main idea of microeconomic the- 
ory, based on independent study of "short-time" and 
"long-time" periods of firm growth. The separation of 
time scales also has deep analogy with methods of study 
of complex physical systems with a wide spectrum of re- 
laxation times, as glasses. For given observation time r 
degrees of freedoms of such systems can be divided into 
"hot" and "cold" ones. Hot degrees of freedoms fluctuate 
in the short-time period (t < t) given that cold degrees 
of freedoms are fixed and can only vary in the long-time 
period (t > r) . Instead of consideration of slow dynamics 
of one system in the long-time period one usually study 
statistical properties of an ensemble of such systems at 
the given time t. 

We apply this approach to find PDF of grow rates 
of firms, which have different dynamics in "short-time" 
and "long-time" periods. In order to establish general 
expression for oscillations of the parameter Ap (t) (10) 
it is instructive to consider first single-harmonic case. 
General expression Ap(t) = V2 a cos (ojt + 4>) can be ex- 
panded over two basis functions £' (i) = cos {tot) and 

(t) =sm(wt): 

Ap (t) - V2a'C' (t) + V2a"i" (t) = V2 (a, £ (i)) . (19) 
which are orthogonal: 

<O = ((0 a ) + ((O a ) = i. <CO = o. (20) 

Here (• • • ) means time average. Instead of two real ba- 
sis functions it is convenient to introduce one complex 
function £ (t) = £' (t) + i£" (t) and complex amplitude 
a = a' + ia" = ae 1 ^, in terms of which the scalar product 
in Eq. (19) is given by expression (a, £) = Re(a*£). In 
the following we use bold notations both for vectors and 
complex numbers. 

In general case, the frequency of quick oscillations 
to > r _1 of £ (t) (as well as its amplitude) randomly varies 
with time. Real and imaginary parts of £ can be con- 
sidered as random values normalized by condition (20), 



where (• • • ) has the meaning annealed averaging over the 
noise £ (t) . Complex amplitude a is fixed in the short- 
time period, and can be considered as random variable 
in the long-time period (or for the ensemble of different 
firms for given time t) . The random function £ (t) and the 
amplitude a describe hot and cold degrees of the freedom 
of the market, respectively. 

Notice, that Ap{t) (19) is invariant with respect to 
"gauge" transformation 



ae 



(21) 



with constant ip, reflecting high degeneracy of market 
quasi-equilibrium states. 



2. Double Gaussian model 



We first calculate PDF of fluctuations Ap, 



V{x)={5 \x-Ap{t)]). 



(22) 



The bar means ensemble (quenched for the time t) aver- 
aging over amplitudes a of fluctuations of different firms. 
The main assumption of "Double Gaussian model" is ex- 
tremely simple: since tactics of firms at the short-time 
period is determined by large number of essentially inde- 
pendent factors, we assume Gaussian statistics of random 
variable £ at time horizon r (due to centeral limit theo- 
rem). But two different firms (or the same firm at two 
different time intervals r) will have, in general, differ- 
ent amplitude of fluctuations a at the long-time strategy 
horizon. Since the strategy of firms is also determined 
by large number of independent random factors, we as- 
sume Gaussian statistics of the random amplitude a with 
dispersion a 2 = a 2 . 

Due to the gauge invariance (21) the noise and the 
amplitude PDFs could depend only on moduli £ = |£| 
and a = |a|. In this section we assume, that hot (£) 
and cold (a) random variables are independent with zero 
average and Gaussian weights 



Qg 



(£) = i e -(«') 2 -(«") 2 , 1 e -[(«') 2 +(«") 2 ]/- 2 (23) 



respectively. 

Fourier transform can be used to calculate the aver- 
ages: 

V(x)= J G (k) e~ lkx ^ G (k) = (ei^k{a'^+a"i"] 

We first calculate the average over Gaussian normalized 
£' and £" and get G (k) = exp{-fc 2 [(a') 2 + (a") 2 }/2}. 
Calculating the average over a' and a", we get G (k) = 
(l + a 2 k 2 /2) . The last step - is to take the inverse 
Fourier transform of this G (k): 
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Exponential distribution of firm grow rates (24) was re- 
ally observed for typical fluctuations x = Ap — —ArG 13 , 
see Eq. (12), with the exponent (3 = 0.15. We conclude, 
that tent-like exponential distribution of firm grow rates 
is the consequence of Gaussian statistics of all degrees of 
freedom (hot and cold) of the market. 



3. Asymmetry of PDF 

The assumption of Double Gaussian model about in- 
dependence of cold and hot variables is, in general, too 
strong, and the noise £ (t) is (anti)correlated with the 
amplitude a. Taking such anticorrelations into account, 
we can write general expression for the noise, satisfying 
gauge transformation (21): 

€(*)=€(<)-Ca/a, a 2 =a^, (25) 

where £ > is the dimensionless correlation factor and 
random variable £ (t) is not correlated with a, and has 

zero average, (t)^ = 0. In the case ( — positive and 

negative fluctuations of firm grow rate, Ar, have equal 
probability, while in the case of positive ( > firms will 
in average grow (because of grow of external resources, 
see section II B). 

At economic level anticorrelations between firm tactics 
and strategy (25) reflect the fact that firms prefer to have 
tactical losses with the hope to get a profit at strategy 
horizons (say, by pressing out business rivals) . And firms 
(and countries), aimed at the maximum instant profit 
without significant investments in the short time period 
will eventually get losses in the long time period. 

Repeating our calculations for the model (25), we again 
find exponential distribution (24) 

„ , . . 1 / e-^/ CT + for x > . . 

Po {x\cr) = ( k (26) 

^(l + C 2 ) I e for x < 

but with different widths a± (a + < cr_) of positive and 
negative PDFs, and the dispersion a: 

°± = " (Vl + C 2 T Cj , ^ = (1 + 2C 2 ) a 2 . (27) 

The average of this distribution is shifted to negative Ap, 
corresponding to systematic tendency to grow: 

JKpj = -V2aC, jA?j = -JKp)G- p > 0. (28) 

Such asymmetrical exponential distribution was really 
observed in the analysis of empirical data in Ref. 27 for 
large averaging intervals (5 years, see Fig. 3). In Fig. 3 the 
x-axis is in units of Ar (12), and not Ap. Empirical value 

1/2 

C = 0.23, and for typical (Ar 2 ) = 0.5 we reproduce 
the observed mean (Ar) = 0.16. 




-2-10 12 
Ar 



FIG. 3: The distribution of grow rates of US firms in 1998- 
2003 for seven size groups from G up = 8—15 through Gdown = 
5 1 2 — 1023 27 . Comparing with the theory we use the same 
correlation factor £ = 0.23 and varied only one parameter 
a ((Ar 2 ) 1/2 = 0.62,0.45 and 0.4,0.3 respectively for upper 
and lower curves. Deviations from exponential dependence 
ar large |Ar| will be explained in section II C 4. 

4- Fat tails 

One of the most prominent features of PDF, the fat 
tail, is usually attributed to large volatility fluctuations 
(in different stochastic volatility and multifractal mod- 
els) . In this section we show, that the tail originates from 
large jumps of the noise, and not of the volatility. This 
new mechanism predicts universal tail exponent fi = 3 
for stock jumps, independent on the coarse graining time 
interval r. 

Fluctuations Ap (t) of the Double Gaussian model are 
characterized by random variable which is Gaussian 
at the time interval r and normalized by the condition 
(£ 2 ) = 1 (20). The problem is that even if wc normalize 
Gaussian variable for given time interval r, this normal- 
ization will be broken at next time intervals because of 
the intermittency effect: relatively rare, but large picks 
of fluctuations. The only way to normalize £ (t) for all 
times is to divide it 

€(*) = €o(*)/Mt). (29) 
by the mean squared average 

(To (t) slowly varies at time interval r, and therefore, ran- 
dom variable £ (t) leaves Gaussian at time scale r. The 
division of £ (t) by ao (t) removes from general Gaussian 
process £ (t) the long-time (at the time scale r) trend 
(long-time variations of the amplitude) , leaving only high 
frequency components. 
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Standard definition of the mean square Op assumes 
that weights Wk in Eq. (30) do not depend on k, and 
we reproduce our previous result (24) for PDF. But this 
definition must be corrected, since there are no any fun- 
damental value of dispersion <Tq, which can only be esti- 
mated from the knowledge of past values of £ . As the 
first step, we have to put Wk = at k ^ and get 

(*) = ™i?o (* - r) + ]T fc>i w k e (t - kr) . (31) 

In the case of totally uncorrelated events ctq is determined 
only by the "reference" value of £q (t — t) at previous 
time interval, and all Wk — ► at k > 1. Terms with 
k > 1 describe the effect of correlations of events, leading 
to variations Ap (i) . 

Second, hot variable £ (t) can vary only on the time 
scale small with respect to r. Therefore, all w k — > at 
k > 2, and random variable £ (t) has Markovian statistics 
with correlations only between neighbour time intervals 
t. Otherwise it will depend on many time intervals time 
kr in the past, which is prohibited by definition of hot 
variable | (t) . 

And the last: the only information known in future 
about past fluctuations, is the very increment Ap, which 
depends only on one component £g = (£oi a ) l a °f Co 
along the vector a. The information about correspond- 
ing "perpendicular" component £ do not enter to the 
increment, and is lost. Therefore, we should drop the 
contribution of (£q) from correlation terms with k > 1 

in Eq. (31): f = + ~^ ■ After a11 these 

corrections we left with expression for the mean square 
in Eq. (29): 

4 (t) = w,e (t-r)+ w 2 % (t - 2t)] 2 (32) 

Although we get similar results for any quickly decay- 
ing weights Wk, calculations are much simplified in the 
case of equal weights w\ — W2 = 1/2 and all Wk — at 
k > 1. In order to calculate PDF of the noise £ (t) (29), 
we rewrite it in the form 

/>oo 

G(0= / da w(a )(6[Z-S /<T }) 
Jo 

f°° a 2 
= / da Tr(a )^e-t 2 °o, (33) 

where we take the average over Gaussian variable The 
probability distribution of the random variable do (32) is 
7T (a) = 2d (6 (cr 2 — Cq)). Using exponential representa- 
tion of this (5-function, we get 

, , a f dse lsa2 [2 2 _sl 

7T (cr) = - / = \ -o l e 2 . 

K J (l + ifi/2) 3/2 V 7T 

Substituting this expression into Eq. (33), we come to 
Student noise distribution: 

Q(£) = ^(l + 2£ 2 r 5/2 . (34) 
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FIG. 4: The distribution of grow rates of US firms in 1998- 
1999 27 , the same parameters as in Fig. 3. Tail exponent /i = 3. 
The varied parameter (Ar 2 ) 1,/2 = 0.45 and 0.3 respectively 
for upper and lower curves. 

Using this distribution function, we finally get 

V(x\a) = -^e^D. 4 (V2-), (35) 
\Jixo \ a J 

where D is the parabolic cylinder function. The central 
part of this distribution has exponential shape (24), while 
its tale has power dependence: 

V (x) ~ |xr 1_M , M»<7. (36) 

with the tail exponent /i = 3, well outside the stable Levy 
range (ji < 2). One can show, that this exponent does 
not depend on relation between weights W\ and u>2 ~ 1 
in Eq. (32) for Markovian noise. But in the absence of 
noise correlations, u>2 — > 0, we get the effective exponent 
fi -> 2. 

If we take into account correlations between the noise 
and the amplitude (see Eq. (25) and discussion therein), 
(£ ) = — Ca/a, and after some calculations we get simple 
expression for PDF: 

/>oo 

V{x) = l da ir (fro) Vo (x\a/ao) = 
Jo 

1 / o+V (x\a + ) for x > 
a^/TTc 1 I aJP fOT X < ' 

where functions Vo (x\cr) and V (x\cr) are defined in 
Eqs. (26) and (35), and a± are given in Eq. (27). We 
show in Fig. 4 that Eq. (37) with £ = 0.23 allows to 
explain both the asymmetry and the shape of empirical 
PDF for different size groups. The size dependence of 
both Fig. 3 and Fig. 4 follows Eq. (12) with exponent 
[3 ~ 0.1 and universal Ap. 

Now we study the stability of the exponent /i for differ- 
ent time periods r. The total increment Ap = y/2(a,£) 
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for two joint intervals r is the sum of corresponding incre- 
ments Api — \f2 (aj, £ 4 ) for each of these intervals. Since 
the amplitude a in Eq. (19) slowly varies on the time 
scale t, we take it the same for both intervals, ai = a 2 , 
and so the noise £ is proportional to the sum of noises 
^ i for these intervals. Each of ^ can be represented in 
the form of Eq. (29), and corresponding dispersions o\ 
and <r 2 depend on the same (but shifted over time) se- 
ries of Gaussian variables £ , Eq. (32). Calculating the 
distribution function of the sum £ = ^ + £ 2 , we find 



Q (£) = J da 1 da 2 J j d Sl ds 2 ^-e-i 2 ^ : 



<7l(72 



1 



it 2 (l + iai/2)(l + t« 2 /2) 



(38) 



where 



y/l+iisx +S 2 )/2' 



1 1 1 



(39) 



The tail of the distribution (38) is determined by small 
a, corresponding to large |si| 3> |«2 1 ~ 1- As the result 
we find, that the distribution Q (£) ~ for the time 
interval 2r is characterized by the same exponent fi = 3, 
as each of for the time interval r. The only difference 
is that this asymptotic behavior can be reached at larger 
£, with respect to the distribution function of 

This observation explains why the fat tail in Fig. 3 
for five year period is shifted to higher |Ar|, with re- 
spect to Fig. 4 for one year period data. Experimental 
observation of the stability of the exponent \i = 3 for 
widely different economies, as well as for different time 
periods 28 r, gives strong experimental support of our the- 
ory. The stability originates in nonlinear correlations of 
the noise, see Eq. (29), while linear correlations vanish, 
((^1,^2)) = 0- To demonstrate the importance of such 
correlations, assume, that the noise £j has tail exponent 
/x, and is uncorrelated at neighboring intervals r. Than 
the exponent of £ ~ ^ 1 + £ 2 f° r the interval 2t is equal 
2fi, and not /x, as follows from our model. 

The systematic study of the distribution of annual 
growth rates by industry was performed in Ref. 29 using 
Census U.S. data. It is shown, that all sectors but fi- 
nance can be fitted by exponential distribution (24). We 
checked the data for finance sector, and show that they 
can be well fitted by Eq. (35) with exponent \i = 3. 



D. Main results 

In this section we considered evolution of the market 
as the result of competition of different firms for external 
resources, by analogy with coalcscent regime in physics 
of supersaturated solutions. This analogy allows to find 
informational entropy of the market, and prove the prin- 
ciple of maximum entropy. 



We demonstrate that in coalescent regime for Gibrat 
mechanism of firm growing the distribution of firms over 
their sizes follows the Pareto power low with the expo- 
nent 7=1 (Zipf distribution). Taking into account size 
effects, it turns to stretched exponent distribution, which 
also describes different processes, related to competition 
of units for common resources. Coalescent mechanism is 
also responsible for observable exponential distribution 
of the income between individuals. The production of 
real firms can be taken into account by vector models, 
by analogy with multicomponent solutions. 

We propose the theory of market fluctuations, based 
on separation of all degrees of freedom of the market into 
cold and hot ones. For Gaussian statistics of all degrees of 
freedom such separation leads to experimentally observ- 
able exponential PDF of firm grow rates. We also prove, 
that this distribution has power tail with universal stable 
exponent fi = 3. 

We find analytical expression for PDF, and show, that 
it reproduces observable shape and asymmetry of the 
distribution of firm grow rates, which is related to exist- 
ing anticorrelations between tactics of firms at short-time 
horizon and their strategy at long-time horizon. In next 
section we apply this approach to study price fluctuations 
on financial markets. 



III. FINANCIAL MARKET 

Dynamics of fluctuations is determined by the spec- 
trum of relaxation times of the system. When all times 
are small with respect to the observation time interval 
t, the state of the market at time t + r depends only on 
its state at previous time t, and dynamics is Markovian 
random process. Short-range correlations of price fluc- 
tuations on the market can be studied using stochastic 
volatility models 30 , but in order to describe real mar- 
kets with multi-time dynamics, the model should take 
infinite-range correlations into account 31 , and has "infi- 
nite" number of correction terms. In addition, to take 
empirically observable excess of volatility into account, 
one has to go at the boundary of stability of such mod- 
els. 

The real market has enormous number of (quasi-) equi- 
librium states and extremely wide spectrum of relaxation 
times, by analogy with turbulence 32 and glasses. Mul- 
tifractal properties of time series can be described by 
phcnomcnological Multifractal Random Walk model 33 . 
Although this model well characterizes scaling behavior 
of price fluctuations, it can not capture correlations at 
neighboring time intervals, which determine "conditional 
dynamics of the market" and can be described by the bi- 
variate probability distribution of price increments 34 . 

In previous section II C we show, that the increment 
of the random value P (t) of the time series 



A T P(t) = P(t + T)-P(t) 



(40) 



has the form of scalar product of two-component random 
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vectors - the noise £ (t) and its amplitude a (t) : 

A T P(t) = V2(a(t),£(t)). (41) 

Hot variables £ (t) vary at the scale small with respect 
to t, while characteristic times of cold variables a (t) are 
large with respect to r. The time r plays the role of the 
effective temperature: at minimal trade-by-trade time, 
r ~ Tfc, the price is almost frozen, while in the opposite 
limit r > t it has random walk statistics. In the in- 
termediate time interval < r < r (of many decades) 
the market has "restricted" ergodicity: only hot degrees 
of freedom are exited, while cold degrees of freedom arc 
frozen and determine the amplitude a of price fluctua- 
tions. 

Here we apply this approach to calculate PDF of price 
increments, as well as various conditional distributions 
and their moments. The dependence of parameters of 
these distributions on observation time r will be studied 
later, in section IV. In section III A we introduce hot 
and cold degrees of freedom of the market. Two simpli- 
fied models are formulated and solved in sections IIIB 
and IIIC. "Markovian" model takes short-time corre- 
lations into account and neglects the effect of long-time 
challenges. "Effective market" model captures such ef- 
fects, but neglects any short-time correlations because of 
trader activity. Although both these models capture es- 
sential part of observable phenomenons of price fluctua- 
tions (extremely small linear correlations - the Bachc- 
lier's first law, "dependence-induced volatility smile", 
"compass rose" pattern 35 and so on), they can not de- 
scribe all the variety of such "stylized facts" 40 . 

In section III D we introduce Double Gaussian model, 
that takes all correlation effects into account, and show 
that it allows to explain the behavior of different types of 
stocks 36 . Analytical solution of this model is derived in 
Appendix D. We demonstrate, that this solution repro- 
duces all observable types of "market mill" patterns and 
gives the mysterious z-shaped response of the market for 
all kinds of asymmetry of bivariate PDF, as well as other 
fine characteristics of this distribution. We also show 
that our theory allows to explain empirically observable 
Markovian "double dynamics" of signs of returns on the 
market 37 . 



A. Cold and hot degrees of freedom 

The idea of hot and cold degrees of freedom of the mar- 
ket is qualitatively supported by empirical observations: 
It is shown in Rcf. 38 , that the amplitude of fluctuations 
for ensemble (quenched) averaging significantly exceeds 
the amplitude of fluctuations for time (annealed) aver- 
aging. This observation can be interpreted as the result 
of the presence of cold degrees of freedom, which remain 
"frozen" when considering time fluctuations of hot de- 
grees of freedom. In the case of ensemble averaging such 
cold degrees of freedom become "unfrozen", increasing 



the amplitude of price fluctuations with respect to its 
time average value. 

Following Rcf. 39 consider two consecutive price incre- 
ments, x (push) and y (response) for the time intervals 
t: 

x = A T P(t), y = A T P(t + r) . 

According to Eq. 41 price increments can be written in 
the form of the scalar products: 

x = V2 y = V2(a 2 ,£ 2 ), (42) 

of complex noises ^ 1 = £ (t) , £ 2 = £ (t + t) and com- 
plex amplitudes ai = a (t) , a 2 = a (t + r) . Complex 
random walk £ (t) in the "tactic" space describes "impa- 
tient" agents. Complex random walk a(t) in the "strat- 
egy" space can be thought of as a result of slow variation 
of composition of the population of such agents on the 
market, as well as the activity of "patient" agents. 

Moduli of complex variables £ f and a^ are normalized 
as: 

<tf) = l, ^ = a\ (43) 

a is the dispersion of price fluctuations 

<A T P 2 (i)) = (A T P 2 (t + r))=a 2 . (44) 

Eqs. (42) are invariant with respect to "gauge" transfor- 
mation of noise and amplitude variables, Eq. (21). 

We will characterize correlations of price increments 
by uni- and bivariate PDFs: 



V (x) = {S[x- A T P (t)]) = / dyV (x, y) , (45) 

V (x, y) = (S[x- A T P (t)} S [y - A T P (t + r)]). (46) 

Using exponential representation of (^-function, these ex- 
pressions can be rewritten in the form 

|V^G(fc,0), (47) 

/°° rlk r°° rln 
_J ^_ e -ik X -i P y G{hph (48) 

where G (k,p) is the Fourier component of PDF 

G(k,p) = (e^rfffl+ipA^ft+Tjy (49) 

The variable y may be interpreted as the response on 
initial push x, which is characterized by conditional PDF 

= IJ^p V{x) = J g e -^G( M ), (50) 
The average conditional response is 

/CO 
dyyV (y\x) (51) 
-CO 

= i f dk _ lkx dG(k,p) 

V{X)J 27T ' dp • 
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FIG. 5: PDF of Russian financial market (finam.ru, 2006) for 
r = 5min (<>- EESR, O - LKOH, □ - RTKM, ♦ - SBER, 
■ - SNGS), solid line shows theoretical prediction (35) with 
/i = 3. 



Averaging the Fourier component of PDF (49) over 
fluctuations of Gaussian amplitude ai = &2 and noise £j 

we find G(k,p) = [l + a 2 (k 2 /2 + p 2 /2 + ekp)] ~\ Cal- 
culating the Fourier transformation of this function (48), 
we get the distribution function 



Vt(x,y) 



1 



'■VT 



>{x 2 +y 2 -2exy) 



2 (1- 



, (53) 



where Kq is the Bessel function. Calculating the inte- 
gral (51) with function (53), we find the conditional re- 
sponse 



(v)x = £x - 



(54) 



Linear dependence (54) with e < well agrees with data 
for Russian market, what can be interpreted as indica- 
tion that Russian investors are oriented only on current 
benefits, mostly ignoring opening possibilities at strat- 
egy horizons. Although linear response (54) is typical for 
ACOR group of stocks with e < (according to classifi- 
cation of Ref. 36 ), this model can not describe essentially 
nonlinear response of other groups of stocks. 



The width of the conditional PDF V (y\x) is character- 
ized by the conditional mean-square deviation 

u\ = J dy(y~(y) x ) 2 V(y\x) (52) 

L_ f dk _ tkx 8 2 G(k,p) 
V (x) J 2n dp 2 p=0 ■ 

Large a x correspond to a large variety of the behaviors, 
the "volatility" . The dependence of a x on x reflects the 
volatility clustering: a x should not depend of x if there 
is no volatility clustering. 

The conditional response (51) and PDF (46) depend 
of correlations between noises and their amplitudes 
in two time intervals. Before formulating general model 
(see section HID), that takes all such correlations into 
account, it would be instructive to study some simple 
limits. 



B. Markovian model 

We first consider the case when the amplitude a (t) is 
not correlated with external challenges at strategy hori- 
zons, and ai = a.2 for two neighboring time intervals. 
We also assume that the noise is not correlated with 
the amplitude, but take into account short range cor- 
relations of the noise, ((d,^)) = £• F° r this Markovian 
model we find Eq. 35 for the probability distribution, 
which describes very well Russian financial market for 
t = 5 min, see Fig. 5. For Gaussian noise we find ex- 
ponential PDF (24) of price fluctuations, which is really 
observed for high frequency fluctuations . 



C. Effective market model 

In general, the amplitude a is varied in response to un- 
predictable external challenges. We first study this effect 
in the model of "Effective market", neglecting correla- 
tions between noise £® in two consecutive time intervals 
r, but taking into account random variations of its am- 
plitude a?: 

= (iM) = *W = K^ (55) 

where v is dimensionless correlation parameter, < v < 
1. As in Markovian model we ignore (anti)correlations 
between noise and amplitude. Correlations of the noise 
£® are induced by trader activity, while the change of 
a stock price in the model of Effective market is deter- 
mined only by an external information, which may be 
considered as uncorrelated random process. 

PDF V [x) of this model is proportional to the 
parabolic cylinder function with power tail exponent 
/i = 2 (see Eq. (36)). Non Gaussian character of the 
noise PDF can be ignored when considering the central 
parts of price distributions, when V (x) takes exponential 
form (24). In order to calculate bivariate PDF, we sub- 
stitute equation (42) in (49) and perform the averaging 
over fluctuations of Gaussian variables a°: 

G(fc,p) = (exp[- ( r 2 fc 2 (^) 2 /2- fe5frt 

«rVK§) a /2-^*P (€?,€§)]>• ( j 

The averaging over noise $® is performed with Gaussian 
PDF (23), and the integral over k and p in expression (48) 
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FIG. 6: Profiles of conditional PDF V (y\x) at different x 
and v — 0.95 (a = $0.04) in comparison with empirically 
observable profiles 42 



FIG. 7: Conditional mean-squared deviation as function of 
x = Ap; the result of Gaussian model with v = 0.95 and 
empirically observable D-smile 43 . 



is calculated expanding the function G (k,p) in powers of 



l=Q y 2l Vi{x)Vi (y). 



(57) 



Here V (x) = V (x) is given by Eq. (24), and functions 
Pi (x) are defined by: 



X 



(58) 



PDF Vq (x, y) is symmetrical with respect to indepen- 
dent transformations of its variables, x — » —x,y — > — y, 
and also with respect to time reversal transformation, 
which corresponds to push-response interchange, x < — > 
y. This function is not analytical in origin, and the ge- 
ometry of equiprobability levels can be approximated by 
\ x \ + \y\ X = const, where A ~ 1 near origin and A ~ 2 
far away from it. 

Profiles of conditional distribution (50) are shown for 
different x in Fig. 6. With the rise of the push x the 
response becomes more flat in origin, in good agreement 
with empirical data. Slight deviations between the theory 
and data at large \y\ 3> a are related to non-Gaussian 
character of the noise (leading to power tail), see Eq. 35 
for more details. Calculating integral (52) in the case 
of Gaussian noise, we get the conditional mean-square 
deviation, 



~2 2 
<7_ = <J 



1 



\v 2 (V2\x\/a 



1 



(59) 



This function is plotted in Fig. 7. It demonstrates the 
so called "dependence- induced volatility smile" ("D"- 
smile), well known from empirical data 43 . At small 
|x| < er the standard deviation of the response (59) 
is smaller than the unconditional standard deviation a, 
while at large \x\ > a it is larger. 

The shape of conditional PDF can also be characterizes 
by the kurtosis, proportional to fourth momentum. One 



can show that, in agreement with empirical data, the 
kurtosis of theoretical conditional PDF Vq (y\x) decreases 
with the rise of |x|. We conclude, that Effective market 
model captures main features of the market behavior, but 
it is enable to describe finite response of real stocks. 



D. Double Gaussian model 

In this section we generalize Effective market model 
to take into account short-range correlations at strategy 
horizon because of activity of traders. Such correlations 
lead to the exponent /x = 3 of the power tail of PDF 
(section II C 4) , and relatively weakly affects the central 
part of PDF: for time interval r about several minutes 
it is estimated as about 5% 39 . We assume normal dis- 
tribution of noise fluctuations ^ for the central part of 
PDF, and neglect the effect of noise-amplitude anticorre- 
lations, which is small at short t, and leads to gain/loss 
asymmetry (see section III D 4) . 

For given noise variables £ i we introduce random vari- 
ables $® of Effective market model, which form orthogo- 
nal basis in the space of random functions see Eq. (55). 
Expanding price fluctuations AP (t) and AP (t + r) over 
this basis, we get: 



A T P(f) = ^(a°,£?) +V2(ei,$ 
A T P(t + r) = v^(a§,*jj)+>/2(e 2 ,tf) 



2 J 



(60) 
(61) 



We consider amplitudes a° of Effective market as Gaus- 
sian random variables, Eq. (55). Non-diagonal ampli- 
tudes £j describe the shift of equilibrium on the market 
because of trader activity. Since there are only two in- 
dependent amplitudes, ai and &2, for two time intervals, 
the amplitudes e\ and £2, can be expanded over two di- 
agonal amplitudes aj and a!]: 



(62) 
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In the case £\ = or £2 = this Double Gaussian model 
is reduced to Markovian model (section IIIB), and in 
the case £\ = £2 = - to the model of Effective market 
(section III C). 

PDF of this model is calculated in Appendix D: 

V (x, y) = 7-o (x cos <fi + — ys'm <f>_, y cos <f>_ + x sin + ) , 

(63) 

where V® is PDF of Effective market model, Eq. (57). 
The distribution (63) depends on only four independent 
parameters: the dispersion a, the correlator of the am- 
plitude v (0 < v < 1), and two angles 0_ and <f> + , de- 
pending on starting parameters {Sy} of our model. The 
correlator v describes the "elasticity" of the market to 
external challenges at the strategy horizon. The angles 
(j>_ and </> + control the feedback between trader expecta- 
tions and real price changes at the tactic horizon. Their 
difference, e = <f> + — <f>_, is taken as small parameter of 
our theory, which controls the correlator of neighboring 
price increments 

(A T P(t) A T P(t + r)) = a 2 s. (64) 

Eq. (63) turns to corresponding expression (53) for 
Markovian model in the limit v — ► 1, and reproduces 
Eq. (57) of Effective market model for tj>_ = <f>, = 0. 

The sum (57) goes only over even 21 because of ne- 
glect of noise-amplitude correlations, (tj,) = 0. In gen- 
eral, there are correlations between noise and amplitude, 
described by a factor £ (see section II C 3). Such correla- 
tions (studied in section III D 4) break the symmetry of 
the conditional average (y) x with respect to positive and 
negative x, and are responsible for the so called Leverage 
effect 44 . 

In our theory we have an hierarchy of small param- 
eters, C < |e| < </> < 1- PDF of Double Gaussian 
model with all nonzero £, (\>. e 7^ has no symmetries 
at all. In the case £ = but 0, e 7^ there is a symme- 
try V (x, y) — V {—x, —y), corresponding to rotations on 
the angle it in the plane (x,y). In the case C = £ = 
but (j> 7^ 0, when there are no linear correlations of 
price (64) at adjacent time intervals, PDF (57) remains 
symmetrical only with respect to mixed transformation, 
V (x,y) — V(—y,x), corresponding to rotations by the 
angle tt/2 in the plane (x, y). The change of sign of y in 
the above equation is related to reversion of the time: on 
reversed time scale one can think about losses in future, 
y < 0, as about gains in the "past". This approximate 
push-response invariance was established first time from 
the analysis of empirical data 42 . And finally, in the case 
C = 4> = e = the function V (x, y) acquires the total 
symmetry x — > —x, y — > — y and x < — ► y of Effective 
Market model. 



1. Market MILL, AC OR and COR stocks 

It is convenient to describe the symmetry of PDF with 
respect to the axes y = by antisymmetric component, 




FIG. 8: Two-dimensional projection of log 4 V a (x, y) with re- 
spect to y = axes a) and y — x axes b) for v = 0.95, 4>_ = 8° 
and <f> + = 8.7°. For comparison sake we show corresponding 
empirically observed pictures c) and d) 36 . 



V a (x,y) = [P (x, y)-V (x, -y)} /2. In Fig. 8 a) we plot 
equiprobability levels of positive part of this function, 
(V a (x 1 y) + \V a (x,y)\)/2. 

For small |g| <C <j> the plot demonstrates four-blade 
mill-like pattern (the "market mill" pattern), that was 
observed first time in Ref. 43 , see Fig. 8 c). To analyze 
these pictures it is convenient to divide the push-response 
plane (x, y) into sectors numbered counterclockwise from 
I to VIII. In agreement with empirical data at e > 
the blades in II and IV quadrants of the (x,y) plane 
are thinner than their counterparts, which extend out of 
I and III quadrants. The situation is reversed at e < 
0. With the rise of |ej the market mill pattern becomes 
distorted and only two corresponding blades of the mill 
pattern left well expressed. 

In Figs. 9 a) and b) we show how the market mill pat- 
tern is deformed with variation of e. Varying the angle 
<f>_ for fixed a, v and <p + , we get good qualitative agree- 
ment with observable patterns, shown in Figs. 9 c) and 
d). We conclude that the theory allows to explain all the 
variety of basic patterns for different stocks 36 , and may 
be considered as the basis for their quantitative classifi- 
cation: Fig. 8 with 4>_ ~ 4> + > (e ~ 0) corresponds 
to the mill pattern (MILL), Fig. 9 a) with 4>_ > + 
(e < 0) corresponds to negative autocorrelation (ACOR), 
and Fig. 9 b) with <f>_ < <j> + (e > 0) corresponds to pos- 
itive autocorrelation (COR). Anti-mill pattern (AMILL) 
with 4>_ ~ cf> + < was never observed in Ref. 36 . 

Similar patterns are obtained for symmetry properties 
of the bivariate PDF V (x, y) with respect to different 
axes y = x, x — 0, or y = —x. As example we show 
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FIG. 9: Changing of asymmetry with respect to y = axes 
with parameter e — <f> + — <f>_. We use v = 0.95 and <j> + = 8° 
and vary the angle <j>_ = 14° &),<(>_= 6° b). For comparison 
sake we show typical patterns observed for different stocks 36 
c) and d). 
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FIG. 10: Probability distribution function for the S&P500. 
Daily data from 31/1/1950 to 18/7/2003 46 a) 5 minute incre- 
ments for 1991-1995 47 b) We show the best fit by Eq. (67) 
with 6 = 0.3 a) and 9 = 0.4 6), for comparison we show by 
dotted lines the best fit by Gaussian PDF. 



Calculating the integral over k in Eq. (50) with this func- 
tion G(k,Q) we find one-point PDF 

V2cr (a{ - oq) 

(67) 

As one can see from Fig. 10 this distribution is in good 
agreement with observable PDF of the Standard&Poor 
500 (S&P500) index, that is one of the most widely used 
benchmarks for U.S. equity performance. 



in Fig. 8 b) equiprobability levels of positive part of the 
function T a (x, y) = [V (x, y)-V (y, x)] /2. The blades 
of this market mill are more symmetric than those in 
Fig. 8 a), in agreement with empirical pictures in Figs. 8 
c) and d). 

An attempt to explain market mill patterns for the 
asymmetry with respect to the axis y — was made in 
Ref. 45 , where "hand-made" analytical ansatz for condi- 
tional PDF was proposed. It was explicitly assumed, that 
the response y depends only on push x at previous time, 
and no long-range correlations were taken into account. 
We do not think, that such Markovian model can give 
adequate description of real market with extremely wide 
spectrum of relaxation times. 



2. Univariate PDF 

Now we calculate one-point PDF, Eq. (47), of Double 
Gaussian model. Expression (D2) of Appendix D for 
the Fourier component G (k, 0) can be represented in the 
form 

G (fc, 0) = (1 + a 2 ajk 2 /2y 1 (l + cr 2 a 2 2 k 2 /2y 1 , (65) 
ol\ = cos 9, a 2 = sin# 



with the angle 8 defined by 



3. Conditional response 

Calculating the integral (51), we find the mean condi- 
tional response 



(y) a 



-sign 



(x) V2a 



2ea\a2 — A e x (%) ~ e 2 (x) 
(a\ - a|) 2 ond (x) - a 2 e 2 (x) 
{eotx — A/a{) ei (x) + (ea 2 — A/a 2 ) e 2 (x) 



(af - a\) [aiei (x) - a 2 e 2 (x)] 



(68) 



sin (26») = VI - v 2 sin (20) 



(66) 



where A — a\a 2 \j (a\ — a\) 2 — v 2 , on and (x) are de- 
fined in Eqs. (65) and (66). In Fig. 11 a) we show how 
the dependence (68) of mean conditional response on 
push x depends on the angle . This dependence has 
zigzag structure for MILL group (e ~ 0), it is almost 
monotonic for ACOR group (e < 0), with linear limiting 
dependence (54)), and is essentially nonlinear for COR 
group (e > 0). Similar calculations of conditional mean 
absolute response (Ij/D^ (which shows how the response 
volatility is grow with the amplitude of the given push x) 
show that in the case e = to a good accuracy it is linear 
in the absolute value of the push (\y\) x — cq + c± \x\, 
in good agreement with empirical data . 

To analyze the asymmetry of PDF V [x, y) with re- 
spect to time reversion in the case of MILL group 
(e = 0), it is convenient to introduce the total incre- 
ment of price during the two time intervals, and also 
the difference of these increments z — 2~ 1 / 2 (x + y) and 
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FIG. 11: The dependence of mean conditional responce on 
push x for different angles (j>_. We use a — $0.04 and cj> + — 8, 
= 0.97 and = 12.5° (ACOR)„ v = 0.9 and tf>_ = 9° 
(MILL), and v = 0.8 and 4>_ = 4.5° (COR). Fitting parame- 
ters were taken to reproduce empirical dependences presented 
for some stocks 36 . Conditional responce (z), for v = 0.95, 
<j>, = 8° and 0_ = 7.7° in line with corresponding empirical 
data are shown in Insert. 



z = (y — x). PDF V (z, z) of these random vari- 

ables takes the form of Eq. (63) with the substitution 
<t> — > if = 7r/4 — (j), describing the rotation of the push- 
response plane by the angle 7r/4. Therefore, both PDF 
and all conditional averages are given by above expres- 
sions under the substitution (y) x — > — (z) z and 9 — ► 9' , 
where 



sin (26>') = y/l - v 2 sin (2<p) = y/l - v 2 cos {2(f) . 

Conditional response (z) also has z-shaped structure, 
and it is shown in Insert in Fig. 11 in comparison with 
empirical data. 

Nonvanishing of average responses (y) x and (z) z 
(Fig. 11) allows one to make some "nonlinear" predic- 
tions (68) about future price changes on the market, 
which can not be obtained from the knowledge of only 
linear correlations: the response y in the next time in- 
terval is correlated with initial increment of the price x 
at small \x\ < a, and is anticorrelated with it at large 
\x\ > a. The order of price increments is also important 
for given total increment \f2~z: for small z < zq ~ a the 
average initial variation x is larger than next one, y, and 
the situation is reverted at large z. 



4- Conditional double dynamics 

In this section we discuss the hypothesis of Ref. 37 , that 
the average return is actually the result of composition of 
two independent signals with Markovian statistics: one 
of them positive, and another one negative. It is pro- 
posed to characterize this effect by average daily returns 
(y~) r and {y+) r given that the previous day had a re- 



X 








a: 




I -0.5 




r c 




IS 


1 

a) 




0.01 0.02 

b) 



FIG. 12: Average daily return given that the previous day had 
a return greater than r c (right) and given that the previous 
day had a return smaller than r c (left). Prediction of Gaussian 
model at 8 — 0.2 a) and empirical data 37 b). 



turn greater than r c and smaller than r c : 

(v-) rc = r dx {y) * v {x) i r dxv {x) ' (69) 

J — GO / J — OO 

/oo / roo 

dx{y) x V{x)j jf dxV(x). (70) 

Calculating integrals (69) and (70) for Double Gaussian 
model in MILL case e = 0, we find at r c > 0: 



(y-)r, = - (f+)r, = 



.4 



a\e x (r c ) 
aiei (r c ) - 



- a\e 2 (r c ) 
otie-i (r c ) 



y/2{a\ - a%) 
aiei (r c ) + a 2 e 2 (r c ) 



V2a 



2 

[ei (r c 



e2 (r c )] 



(71) 

This function is shown in Fig. 12 in line with empirical 
data. As one can expect from Fig. 12 a) the average re- 
sponse (y+) r is correlated with r c at small r c < a and is 
anticorrelated with it at larger r c . We added horizontal 
dotted line in Fig. 12 b), shifting the y-axis by uncon- 
ditional average return 37 (y) — 0.00025. This shift and 
remaining difference in shape between Figs. 12 a) and b) 
is related to the buy/sale asymmetry, discussed below. 

At r c = the average daily return for given sign of 
the previous day return is finite, reproducing the effect 
of "double dynamics" of the market, attributed in Ref. 37 
to "propagation" of two independent signals in "Marko- 
vian" market. 

In fact, the market is not Markovian, but the sign 
of price increments is determined only by the noise £. 
Markovian "double dynamics" of signs is direct conse- 
quence of Markovian statistics of noise correlations, see 
section II C 4. Anticorrelations between the noise and the 
amplitude are responsible for small systematic trend of 
the price, (y) ~ \/2aC, (28), reproducing empirical data 37 
for C ~ 0.02. This positive trend leads to corresponding 
increase of the probability to have positive price incre- 
ment, p+ = 1/2 + ciC with ci ~ 1. Conditional probabil- 
ities of the two-state model 37 can be expressed through 
the angle <j) ~ <p_ ~ <f> + , which determines the ampli- 
tude of the response (y) x on previous price increment x 
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FIG. 13: Theoretical dependence of the skewness p x on x for 
double Gaussian model at <f) = 0.2, e = and v — 0.9 and 
empirically observable dependence 39 . 

(ca ~ 1): 

p ++ = 1/2 + ciC + c 2 0, 
P— = l/2- Cl C + c 2 0. 

Empirical observation of "double dynamics" may be 
considered as direct confirmation of Markovian statistics 
of noise fluctuations, but not of the whole market, as 
conjectured in Ref. 37 . We show in section IV C 7, that 
multifractal evolution of the amplitude is strongly non- 
Markovian. But consideration of only signs of returns 
"erases" information about the amplitude from the time 
scries. 



5. Skewness 

Asymmetry of the conditional distribution V (y\x) 
with respect to the average (51) is characterized by the 
skewness of the conditional response: 

1 f°° 

Px = -3 / d y (y - (y) x ) ? (vV) ■ 

The conditional mean-square deviation a x is defined in 
Eq. (52). Positive value of p x indicates that only few 
agents perform great profits, while many of them have 
small losses with respect to the mean. A negative p x 
describes a complementary case. As one can see from 
Fig. 13, the skewness has the sign of initial push x in 
accordance with the empirical dependence. Notice that 
although the skewness is very sensitive characteristic of 
PDF, our theory reproduces both observed shape and 
values of p x . 

In this section we show, that separation of hot and 
cold degrees of freedoms allows to reproduce numerous 
empirical data, known for high-frequency fluctuations on 



financial markets. For Gaussian statistics of all degrees of 
freedom this model captures main features of all groups 
of stocks, including "market mill" patterns, "dependence- 
induced volatility smile" , z-shaped response function and 
so on. Correlations between hot and cold variables are re- 
sponsible for observable double dynamics of the market, 
mixing propagating signals of opposite signs and provid- 
ing systematic positive trend of prices. 



E. Results and restrictions 

In this section we demonstrate, that the idea of hot and 
cold variables allows to capture main features of price 
fluctuations, which can be described by Double Gaus- 
sian model - a generalization of the random walk model 
for the case of multiscale fluctuations. For different sets 
of parameters the analytical solution of this model re- 
produces the behavior of all kinds of stocks on financial 
market, as well as the market as a whole. 

Consider some restrictions of this approach. Using 
coarse-grained description of price fluctuations at times 
T > T k ^ 1 noin we loose information about sale/buy 
mechanisms 48,49 . This knowledge (see, for example, Mi- 
nority and Majority Games 50 ) is important to derive pa- 
rameters of our models for particular markets. We also 
considered only uni- and bivariate distribution functions, 
while market dynamics is described by the whole fam- 
ily of rt-point correlation functions. In next section we 
present alternative description of the market at differ- 
ent time scales using ideas of renormalization group ap- 
proach. 



IV. MULTISCALE DYNAMICS OF THE 
MARKET 

Standard thermodynamics can describe only ergodic 
systems, while the market is the system with "restricted 
ergodicity" : during the time r it can explore only small 
part of the total configuration space near current local 
equilibrium. Increasing the time t, this equilibrium is 
shifted because of long-time variations. The resulting 
multi-time dynamics of fluctuations on the market is not 
ergodic, and can only be described by continuous set of 
Langevin or Fokkcr-Planck-type equations for all time 
scales. 

Note, that this is not exclusive, but standard behavior 
of complex physical systems. As we will show, different 
local equilibriums on the market are organized into tree 
cascades ("self-organized criticality" ) , by analogy with 
"hot spots" in Quantum Chromodynamics 51 , and dy- 
namics of unergodic spin-glasses, which is government 
by continuous set of Fokker-Planck equations 52 . Contin- 
uous set of equilibriums was also predicted for incomplete 
markets 53 . 

The behavior of the market at the trade by trade level 
was studied in many details 54-56 . At larger times col- 
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lective effects become important, and financial time se- 
ries display long-time nonlinear correlations 57-59 , which 
puzzle many researchers 60-62 . Different models have 
been proposed in order to reproduce some "stylized 
facts" 40 of empirical time series. Levy flight processes 63 
were used to model jumping character of price varia- 
tions. Volatility (the amplitude of fluctuations) cluster- 
ing effects have been studied in frameworks of stochastic 
volatility models 64 and GARCH-type models 65 . Mixed 
effect of jumps and stochastic volatility was taken into ac- 
count in some models 66 . A key to study multiscale prop- 
erties of price fluctuations is provided by phenomcno- 
logical multifractal models, see Refs. 67,68 . Renormgroup 
approach, describing evolution of multiscale systems, was 
first proposed for glass systems 69 ' 70 , and in present work 
we extend it to the market. 

The key question is the source of price fluctuations. 
Assumption about totally random activity of traders 71 ' 72 
leads to Brownian motion of prices. Although random 
trading model predicts many qualitative and quantita- 
tive properties of the order books 55 ' 73 , it can not describe 
existing correlations in price fluctuations. An alternative 
"efficient market hypothesis" assumes, that the price can 
be changed only because of unanticipated and totally un- 
predictable news. This hypothesis lays in the basis of 
the model of fully rational agents 74 , which also predicts 
Brownian walk statistics of prices. Observed volatility of 
the market is too high to be compatible with the idea of 
fully rational pricing 75 , and can only be reproduced by 
introducing an artificial random source - "sunspots" 76 . 
In addition, the analysis of Ref. 77 shows, that most of 
large fluctuations in the market are due to trading activ- 
ity, independently of real news. 

The main idea of this paper is that market activity can 
be described as random trading at all time horizons t p 
from a minute to tenths years. The market tends to reach 
equilibriums on extremely wide baseband of time scales, 
and all these equilibriums are continuously changed both 
because of long-time modes and external events. Multiple 
local equilibriums can be represented by an hierarchical 
tree, see Fig. 14. Each generation r of this tree is char- 
acterized by its relaxation time r r . For any observation 
time interval r — r r , all states with times smaller than 
r r are in equilibrium, and fluctuations near this equi- 
librium are described by hot degrees of freedoms. The 
states with times larger t t are frozen, and are described 
by cold degrees of freedom. 

In this section we derive an analog of renormgroup 
equations for the market, which relate fluctuations at 
different coarse grained time scales r. With decrease 
of the time interval r, which can be thought of as ef- 
fective temperature, the market experiences a cascade 
of dynamic phase transitions of broken ergodicity, when 
some hot degrees of freedom become frozen. This cascade 
can be graphically shown as the hierarchical tree, each 
branching point of which represents "phase transition" 
to a state with frozen degrees of freedom with relaxation 
times T r > t, see Fig. 14. We show in section IV A 1, 
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FIG. 14: Hot and cold degrees of freedom have, respectively, 
times small and large with respect to coarse graining time r. 
The scale of relaxation times for two observation times r and 
t/ (/ — 1), Fig. a). Elementary step of renorm group trans- 
formation corresponds to division of "parent" time interval 
t into (/ — 1) "child" time intervals rj (/ — 1), see Fig. b). 
Hierarhical tree in ultrametric time "space" is shown in Fig. 
c). For given observation time r upper part of the tree, shown 
by solid lines, corresponds to "cold" degrees of freedom, and 
lower (dotted) part corresponds to "hot" degrees of freedom. 



that the topology of this tree reflects ultrametricity of 
the time "space" . 

In section IV A 2 we demonstrate, that fluctuations at 
given time scale r are determined by contributions of all 
"parent" time scales of the hierarchical tree in Fig. 14, 
what is the reason of non-Markovian dynamics of the 
market. Cumulative contribution of all time scales al- 
lows to explain extremely high volatility of the market 
(section IVB1), and is responsible for power low de- 
cay of correlation functions (section IV CI) and their 
multifractal properties. In section IV C 1 we formulate 
self-consistency condition, under which the hierarchical 
tree in Fig. 14 describes coarse-graining dynamics at all 
levels of the coarse graining time r, and find the r- 
dependence of parameters of our Double Gaussian model, 
section HID. 

In section IV C 3 we derive a set of Langevin equa- 
tions, describing dynamics of the market with extremely 
wide range of characteristic times - from minute to tenths 
years, and calculate the price shift in the response on im- 
balance of trading volumes (section IV C 4). 

We also calculate PDF of volatility (section IV C 5) 
and show, that it has fat tail with stable exponent /j = 3 
for stock jumps and /j = 2 for news jumps. We de- 
rive, that coarse grained dynamics of the market can 
be reduced to the multifractal random walk model 78,79 , 
which determines multifractal properties of price fluctu- 
ations, related to the ultrametric structure of the tree 
in Fig. 14. We calculate volatility patters after news 
and stock jumps, and find their conditional probabili- 
ties. In section IV C 6 we show, that the price P (t) be- 
haves as fractional Brownian motion. We demonstrate 
in section IV C 7, that Brownian motion, sub- and super- 
diffisive regimes change each other at the long-time scale. 
The knowledge of history can be used to estimate the 
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tendency and risks of future price variations. 

Main results of this section are summarized in sec- 
tion IV D. In Appendix E we show details of calculations 
of the volatility distribution. 

A. Renormalization group transformation 

Consider statistics of price increments (returns) 

A T P(t)=P(t + r)-P(t), (72) 

as the function of the coarse-graining time interval r. 
Here P (t) is the price or its logarithm at time t. For 
definiteness sake, we consider only ACOR group of 
stocks 36 , when A T P (t) can be represented as scalar prod- 
uct of complex amplitude a (t) and complex noise £ (t) , 
Eq. (41). By analogy with renormgroup consideration, 
cold variable a (t) slowly varies at time scale r, while hot 
variable £ (t) quickly fluctuates at this scale, see Fig. 14 
a). 

1. Ultrametricity and restricted ergodicity 

In order to establish an analog of renormgroup trans- 
formation for the market we first introduce correspond- 
ing partitioning of the time "space". At elementary 
step of the renormgroup each "parent" interval r r of 
time axis can be divided into / — 1 "child" time interval 
T r +i = T r / (/ — 1), see Fig. 14 b). Repeating such divi- 
sion, we arrive to the hierarchical tree with functionality 
/, shown schematically in Fig. 14 c). For this tree the 
time 

T r = T e~ Kr , K = ln(/-1) (73) 

depends exponentially on the current rank r, To is maxi- 
mal relaxation time. Minimal time = Toe~ Kk (k is the 
number of generations of the tree) is about average time 
between trades. Typically, To is about several years and 
Tk is about minute for liquid markets, and so nk > 13. 

We define the "distance" z between events at times t 
and t' by the condition t,._ z = \t — t'\: 

z (t - t 1 ) = - In ^ ~ 1 ^ at \t-t'\ >r = r r , (74) 

which can be identified with the distance (number of gen- 
erations) along the tree in Fig. 14 c) between these points. 
One can show, that for three different times t,t',t" 

z(t- t") ~ max [z (t -t'),z (t' - t")} , 

and, therefore, the metric (74) generates ultrametric time 
"space" (with only isosceles and equilateral triangles), 
which can be really mapped to the tree. For example, in 
Fig. 14 c) the distance between points a and b is z ab = 1, 

z bc = 2, and z ac = max (z ab , z bc ) = 2. 



Note, that Eq. (73) gives standard relation between 
time scales of discrete wavelet transformation. The tree 
in Fig. 14 can be thought of as a skeleton of the wavelet 
transformation of time series. We turn to wavelet inter- 
pretation of our approach in section IV C 3. 

Each of horizontal levels at the time r = r r on the tree 
in Fig. 14 c) corresponds to course-grained description of 
fluctuations at the time scale r. Hot degrees of freedom 
are "melted" , and described by complex noise £ (t) with 
continuum spectrum of relaxation times extended from 
Tfc through t. Cold degrees of freedom are characterized 
by complex amplitude a = a r of the noise, which is frozen 
at the time r, see Eq. (41). 

By analogy with glasses, the states of real market 
are highly degenerated, what is reflected in the pres- 
ence of gauge transformation (21) of complex noise and 
amplitude variables, which do not affect price variations 
A T P (t), Eq. (41). The degeneracy is the reason of high 
sensibility of the market to external events. Recall, that 
in spin glasses any observable are not affected by "non- 
serious" part of disorder, which can be removed by gauge 
transformation of glass degrees of freedoms. 

Following this analogy, the time r plays the role of the 
temperature T. With decrease of the temperature T ~ r 
from r = To the market experiences a cascade of dynamic 
phase transitions of broken ergodicity, when some hot de- 
grees of freedom become frozen (the system is unergodic 
if its fluctuations can not explore the whole configura- 
tion space) . This cascade proceeds continuously down to 
the time Tfc, and can be graphically shown as hierarchi- 
cal tree, each branching point of which represents phase 
transition to a state with frozen degrees of freedom with 
relaxation times r r > t, see Fig. 14. The parameter 
k <C 1 determines the probabilities of such transitions, 
which are relatively rare for real markets. 

In this sense at any t < To the market is just at the 
point of dynamic phase transition of broken ergodicity, 
and has, therefore, increased amplitude of fluctuations - 
the volatility. This observation supports the idea that the 
market is always operating at a critical point as the result 
of competition between two populations of traders: "liq- 
uidity providers", and "liquidity takers" 80,81 . Liquidity 
providers correspond to hot degrees of freedom, creating 
antipersistence in price changes, whereas liquidity takers 
correspond to cold degrees of freedom, and they lead to 
the long range persistence in prices. 

Since such separation of the market into hot and 
cold degrees of freedom takes place at any time scale, 
Tfc « t « To, there could not be any unique classifica- 
tion of traders, which can be divided also into "positive 
feedback" traders and "fundamentalists" 82 , "contrarian" 
traders and "trend followers" 83 and so on. There is, how- 
ever, important difference between market and spin-glass 
hierarchical trees: while the states of the glass are not or- 
dered, there is strong time ordering of all "points" of the 
market tree at any level r of the hierarchy. 
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2. Recurrence relation 

General recurrence relation between amplitudes a r and 
a r+1 at levels r and r + 1 of the hierarchical tree can be 
written through the random transition matrix u r : 

a,-+i = u r+ ia r + Aa r+ i. (75) 

In general, there could be a term ~ (a r )* in the rhs of 
this equation, but it is not invariant with respect to gauge 
transformation, Eq. (21), and should be dropped. The 
term ~ u r+ i describes the inheritance of the amplitude 
a r of the "parent" levels of the hierarchy, and Aa r+ i 
gives the contribution of "newborn" during the transition 
r — > r + 1 unfrozen degrees of freedom to the "child" 
amplitude a r+ i. 

The recurrence relation (75) can also be rewritten in 
the multiplicative form, introducing relative increment 
A r+i = Aa, r+1 /a r : 

a r+ i = e" r+1 a r , e" r+1 ~ l+w r+ i = u r+ i + A r+ i (76) 

Random variables u r+ i and Aa r+ i (A r+ i) are deter- 
mined by degrees of freedom with characteristic times 
r r < t < r r +i, while a r is formed by degrees of freedom 
with times larger r r . Assuming independence of fluctua- 
tions of different time scales, a r do not depend on u r+ i 
and Aa r+ i. We estimate the mean squared amplitudes 
of fluctuations of Aa r and u r as 

(Aa r ) 2 = D T r , tf r = u 2 , (77) 

There is important difference of Eq. (77) from the case 
of usual diffusion, when (r 2 ^ = Dt is the consequence 
of independence of fluctuations at different times t. In 
contrast, diffusion-like dependence (77) with effective co- 
efficient D is the consequence of independence of fluctu- 
ations of different time scales r r . The time t-dependence 
of price fluctuations is strongly non-diffusive. 

B. Amplitude of fluctuations 

1. Excess of volatility 

In the mean field approximation we neglect fluctua- 
tions of u r = it, and find the solution of Eq. (75) in the 
form of the sum of independent random signals Aa r _fe 
from time intervals r (/ — l) k , obtained by multiplicative 
merging of (/ — 1) previous time intervals r: 

^ W = E fe>0 ukAa r-k (t) . (78) 

Weights of these signals exponentially fall with the dis- 
tance k in time hierarchy from the current rank r. Simu- 
lated time series (78) for the amplitude a (t) in the model 
with random Aa r = ±\/D r r are shown in Insert in 
Fig. 15. This picture demonstrates multiscale character 
of resulting price fluctuations. 
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FIG. 15: Empirical dependence 84 of dispersion a (r) on time 
interval r, and its fitting by Eq. (82), Aq = 0.9 for DELL and 
Aq = 0.8 for General Electric. The effective Hurst exponent 
H is different at r < t x and r > t x . In Insert - computer 
simulation of random amplitude, Eq. (78). 

Averaging the square of the recurrence relation (75), 
we find difference equation 

a 2 (r r+1 ) = u 2 a 2 (r r ) + D Q T e~ K ^ (79) 

for the dispersion a (i>+i) of the amplitude a r , which 
has the solution 

a 2 (j r ) = Dr r + Lu 2r , (80) 

where L is the constant of integration and D is the effec- 
tive diffusion coefficient: 

D = ^— - > D , e- KX ° = u 2 e K . (81) 

1 — e o 

From Eqs. (80) and (73) we find the dependence of 
the dispersion of price increments on the coarse-graining 
time r: 

ct 2 {t) = Dt + L(t/t ) 1+x1 . (82) 

The dependence (82) for different stocks is in good agree- 
ment with empirical data, see Fig. 15. Although at small 

T < T x , 

t x =t (Dt /L) 1/X \ (83) 

it looks like diffusion with apparent diffusion coefficient 
D (81), price fluctuations do not really have diffusive 
behavior. As the sign of it, the amplitude of price 
fluctuations is anomalously large due to the presence 
of a big prefactor for k < 1 in Eq. (81). It was 
shown by Schiller 85 , that even accounting the volatility of 
dividends 86 leaves the empirical volatility at least a factor 
5 too large with respect to the random walk model. Such 
anomalous "excess of volatility", D/Dq ~ 10, originates 
from the superposition of signals from all time scales, see 
Eq. (78). Similar effect (by 10 orders of value) is well 
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known in spin glasses, where the parameter k in Eq. (74) 
is extremely small. 

Eq. (80) can be used to estimate the amplitude of the 
transition matrix in Eq. (76). From Eqs. (77) and (80) 

with L = (at t < t x ) we find A 2 +1 = er K (l - e~ KA o) , 

and get 

7^=e- K . (84) 



2. Cross-over time and Hurst exponents 

The term ~ L in Eq. (82) appears as a constant of 
integration of the recurrence equation (79), which is de- 
termined by the "boundary condition" at trading time 
Tfc. Therefore, L is not universal and determined by mi- 
crostructure of the market. 

At small t < t x the first term in Eq. (82) gives the 
main contribution, a (r) ~ t h , with effective Hurst ex- 
ponent close to 1/2. At large r > t x the second term 
in Eq. (82) gives dominating contribution to the Hurst 
exponent: 

H=(l + Ag)/2. (85) 

Such behavior with different exponents H at r < t x 
and t > t x was really observed for S&P 500 stock in- 
dex (1984-1996) 87 with different values of the cross-over 
time t x for individual companies, see Fig. 15. 

The removal 87 of the largest 5 and 10% events kills 
correlations of the noise £ (t) at small time scales, reduc- 
ing the constant L. According to the prediction (83) of 
our theory, it shifts a (r) to lower values, and strongly 
increases t x . Excluding the shift of L from variations of 
both a and t x , we find linear relation between two these 
shifts: Alncr ~ — (Aq/4) Alnr x . Comparison with em- 
pirical data 87 gives the estimation Aq ~ 1 , in good agree- 
ment with observable exponent H ~ 0.93 for the regime 
t > t x , see Eq. (85). Similar behavior is observed for 
different stocks with typical transition times t x about 
several days. 



3. Parameters of Double Gaussian model 

Let us show, that in order to represent coarse grained 
dynamics of price fluctuations for the time interval t = 
T r , noise variables at different "child" subintervals 
n = 1, ...,/— 1 of the same "parent" interval should 
be (anti)correlated. According to the idea of the coarse 
grained description, the price increments for the time in- 
terval r r is the sum 

A Tr P(t) = J2 f ~ 1 A Tr+1 P(t + nT r+1 ) (86) 

of price increments for all / — 1 adjacent time intervals 
7>+i. Substituting Eq. (41) into Eq. (86), this last equa- 



tion can be rewritten in the form 

( a r,£r) = £ ( a r,£r+l) 

(/-I \" /-I 

u r a r ,^C+i +£( Aa ?+i'£n)> 
n=l / n=l 

where we substituted Eq. (75) for the amplitude a" +1 to 
the right hand side of Eq. (87). Calculating the aver- 
age (both quenched and annealed) of the square of this 
equation, we find the self-consistency equation 

a 2 ( Tr ) = U V (T r ) (f - 1) (1 + s) + D T r (f - 1) , 

where e is the noise correlator at neighboring time inter- 
vals, 

£ = (A Tr+1 P(i) A Tr+1 P (t + T r +l)) j <J 2 ( Tr+ l) 

= <(e 1 ,^ +1 )>' 

Substituting here u 2 = e _K ( 1+A o) (80) with e K = f — 1 
we find in the case L = 

e ~ - (e K - 1) k\ 2 q ~ — k 2 Aq. 

This value is always negative and small at n <C 1. 

4- Time and size dependence of fluctuations 

The effect of noise/amplitude anticorrelations, studied 
in section II C 4, is small in the parameter ( C 1, and 
leads to the asymmetry of the tails of probability dis- 
tributions (see Fig. 4), observed for PDF of individual 
companies in Rcf. 88 . It is also responsible for different 
apparent exponents for positive (/i + > 3) and negative 
(ji_ < 3) power tails. This effect is illustrated in Insert 
in Fig. 16, where we show, that the function (x — x) 4 
at x = 1 is indistinguishable at about two decades in x 
from power tails \x\ 1 M± with exponents /i + = 3.25 and 
fi_ = 2.8. Increasing x ~ C leads to larger deviations of 
these exponents from the universal value 3, see Insert. 

Simple analytical expression for the whole distribution 
V (x) for general r can be easily constructed by consec- 
utive matching Gaussian, exponential and power distri- 
butions, V (x) ~ x -4 , at some points x± and y±, respec- 
tively. The resulting expression well reproduces most of 
empirical data for PDF V (x) . Below we use this expres- 
sion to explain main features of PDF at r <§C t x and 
t > t x : 

We can estimate the effective exponents fi ± re- 
expanding x~ 4 ~ (x — x) 1 M± about correlation in- 
duced systematic shift x ~ ^/2C,a ~ V2(cr (see Eq. 28), 
and find [i± ~ 3 + 4x/x± near cross-over points x = x±. 
At small r <C t x the central part of the distribution has 
exponential shape (26) with maximum at x. Matching it 
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FIG. 16: Empirically observed positive a) and negative b) 
tails 88 of the probability distribution, and their matching by 
shifted Gaussian and power tail (x — x)~ 4 , x = O.ler for time 
interval r = 64 days and x = l.lcr for r = 1024 days. In Insert 
we show matching of x~ 4 by (x — x)~ 1 ~'~ l± with = 3.25 
= 2.8 (2) at x = 1, and ^ + = 4.5 (3),/Lt_ = 2.3 (4) at 

x = 2. 



at a; = a;± with power tails A±x 4 we find x± — ±2\/2cr± 
and fi ± ~ 3 ± 2£, in good agreement with Ref. 88 . 

At large r S> t x the central part has Gaussian shape. 
Matching it at x = y± with power tails, we get y± = 
x/2± y/x 1 /4 + 4cr 2 . With increase of C (at x ~ \/2Ccr > 
ct) Gaussian region of the positive distribution is pro- 
gressively extended, while negative distribution remains 
fat-tailed, explaining corresponding mysterious behavior 
of empirical data 88 , see Fig. 16. 

The size dependence of the volatility was studied for 
individual companies in Ref. 88 . It is shown, that Eq. (13), 
<7 ~ G - ' 3 , well describes the dependence of dispersion of 
returns on market capitalization. For r = 1 day (3 ~ 0.2, 
while it gradually decreases with the rise of t, approach- 
ing the value (3 ~ 0.09 for r = 1000 days. This effect 
supports our self-similar model of companies (see sec- 
tion II B 2), when the index (3—1/ (2n) is determined 
by the number n of generations of the hierarchical tree. 
The effective number n of tree generations logarithmi- 
cally depends on relaxation time r, Eq. (73), and for 
Tfc <C t <C To the dependence (3 on r can be approxi- 
mated by 



/^Aj-Alnr. 



From equation er ~ G ^ ~ r ff we find that the Hurst 
exponent (85) should grow logarithmically, 

= H + /3 1 mG, 



with market capitalization G, in good agreement with 
empirical data 84 . 



C. Nonlinear dynamics of fluctuations 

1. Correlation functions: multifractality 

From Eqs. (75) and (80) we find simple analytical ex- 
pression for the correlation function of amplitudes: 

(a (t),a(*')) = £>r e-' cr -' a ° z ( t > t ') + Lu 2r , (89) 

where z is the logarithmic distance (74) in the ultrametric 
space, see Fig. 14. Therefore, observed power autocorre- 
lations in the time series are the consequence of the self- 
similiarity of the hierarchical tree in Fig. 14. Neglecting 
the term L at large enough r (small coarse-graining time 
r < t x ) in Eq. (89) we find that amplitude correlation 
function decays as the power of the time 



(a (t) , a (tf)) ~ exp (-A 2 , In \t-t'\ /r) , \t - t'\ > r. 

Now consider fluctuations of the modulus a (f ) of the 
amplitude a(t). The solution of the multiplicative re- 
current relation (76) for the coarse graining time t = r r 
is 



a (i) ~ doe' 



;(t) 



(*) = E <r W P (*) ' 



(90) 



From expansion (90) we find 



a? (t) a? (i') 



r(9) , T(?) = ^M_^M 



(91) 



with 



■ (?) = Lae^- = -g/s/2 + kA 2 (q 2 /2 - q) + 



where we expanded g (q) over irreducible correlators of 
uj r = uj r + Aw r , kA 2 = Aw^, and used Eq. (84) to find 
the linear in q term. For Gaussian uj r there are only two 
first terms in this expansion, and we get t (q) = X 2 q 2 . 
We also find 



a i (t) ~ T qH W , gff (g) = -g (g) /«. 

In the case of Gaussian w p this gives us the generalized 
Hurst exponent H (q) = 1/2 + A 2 - X 2 q/2, see also Ref. 89 . 
The intermittence parameter A 2 characterizes the uncer- 
tainty on the market, and we expect, that A 2 is relatively 
large for emerging markets with large uncertainty, and 
small for well-developed markets (see section IV C 3). 

We conclude, that hierarchical structure of market 
times, see Fig. 14, generates multifractal time series with 
g-dependent generalized Hurst exponent. The ampli- 
tude a of fluctuations is randomly renewed with time 
t: with the probability p ~ Tq 1 for the root of the hi- 
erarchical tree in Fig. 14, and with the probability 
Pk ~ 1 S> 7,. for maximum rank i = k of the tree. 
This random process generalizes the Markov-Switching 
Multi-Fractal process 67 with a 2 = cr 2 ]J k r=1 M^. The 
multiplier is renewed with probability p r exponen- 

tially depending on its rank r within the hierarchy of 
multipliers. 
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2. Volume statistics 

In this section wc introduce an analog of canonical 
action-angle variables, in which coarse-grained market 
dynamics can be described by a set of linear Langeven 
equations for all time scales t p in the market. The "ther- 
modynamic force" of price variations is the imbalance of 
trading volumes, V (t), which can be considered as ran- 
dom function of time (volume time series). The incre- 
ment of the volume at the time interval r = r r can be 
written by analogy with price increment (Eq. 90) in the 
form: 



A T V(t)~a e v Wr)(t), v(t) 



v p (t) . (92) 



Normalized random noise rj (t) is proportional to the sign 
of the increment A T V (t). Gaussian random variable v (t) 
slowly varies at the time scale r, and can be expanded 
over modes p covering the frequency band from r~ 1 to 
Tp+i- Explicit expression for v p (t) can be obtained ex- 
panding its variation Av (?) = v p (?) — v p over normalized 
wavelet functions ip with expansion coefficients v p (?) (v p 
describes regular trend): 



Av p (t) = Jt p 1/2 i, [(t - if) /t p ] v p (?) dt', 

% (?) = I Tp 3/2 r [(t - t') /Tp] AV (t) dt. 



(93) 



Similar equations relate modes uj p (t) (90) with the 
volatility oj (?) . In the case of random activity of traders 
v p (if) at all time horizons t p should be considered as 
uncorrelated random values: 



v p (t) v p i (?) ~ A k6 pp i5 (t — ?) . 



(94) 



The noise function r\ (t) in Eq. (92) can be presented 
in the form rj (t) ~ cos (t) , where <j> (t) is the phase 
of corresponding complex amplitude (see Eq. 41). Ran- 
dom function tfi (t) can be expanded over wavelet modes 
(Eq. 93). Assuming that corresponding expansion coeffi- 
cients <j) p (?) are independent random values at all time 
horizons, 



cj> p (t)<P p , (?)^ 1 K8p pl 5(t-?) 



(95) 



wc find: 



r] (?) r) (?) - r] (?) ~exp 

J Tk kt p 



E <a p (*)^(*') 

k <T p <\t-t'\ 



~ exp 



Calculating the integral over t p we arrive to 



r)(t)r)(?)-r,(t) 



t-? 



(96) 



From above equations we find that the amplitude of 
volume increments | A T V\ is log- normally distributed and 



uncorrelated with the sign of A T V (t). These predictions 
are supported by empirical data 90 , which also show, that 
signs -q (?) of trade volumes (and, therefore, the very -q (t)) 
have long-range power correlations, Eq. (96), with stock 
dependent exponent 7 < 1. We conclude, that observed 
power-low correlations in signs of volume are the con- 
sequence of the self-similiarity of price fluctuations at 
different time scales, which lead to scale invariant in- 
termittence parameter A 2 (Eq. 94) and the amplitude of 
fluctuations 7 (Eq. 95). Such long range correlations are 
usually considered as the result of cutting of large trades 
into small chunks (see Ref. 90 ). 



3. Langeven equations and market entropy 

The key idea of Langeven formulation of multi-time 
market dynamics is that fluctuations at different time 
scales t p are statistically independent. Therefore, the 
logarithmic volatility io p (t) of the mode with relaxation 
time t p is induced only by corresponding volume mode 
v p (t) (93). Since both uj p (t) and v p (t) have Gaussian 
statistics, general Langeven equations are linear (differ- 
ent choice of fluctuation modes makes these equations 
highly nonlinear): 



du> p (t) 



p dt 

with (5-correlated noise: 



+ u p (t) = v p (t) 



Av p (t) Av p/ (?) = 25pp,n\ l T p 5 (t - ?) . 



(97) 



(98) 



At time scale t p this equation is in agreement with 
Eq. (94). After standard calculations we find correlation 
function of volatility modes 



Auj p (t) Aup, (?) = 8pp,K\ 2 e-\ t -~ t '\ /T e, 



(99) 



and fluctuations of logarithmic volatility for the coarse 
graining time r = r r : 



G(t — ?) = Aoj (?) Aoj (?) - £ p ^ r Aoj p (t) Aoj p (?) 
~ A 2 ln(T /|i-i'|), t «f|i-i'|«T . 

( 10 °) 

This expression lays in the basis of multifractal random 
walk model 78 ' 79 , and reproduces Eq. (91): 

ai (t) a* (?) - exp<z[w (t) + oj (?)} - (r/ \t - t'\) T(q) . 

Eqs. (90), (97) and (98) present Langeven formulation 
of multifractal market dynamics. Using standard trans- 
formations, they can be rewritten in the form of Smolu- 
chovski equations for probability function ^ {w p }. The 
importance of this function is that it defines rigor entropy 
of the market 

S[*] = J Dw p ^> {w p }ln* {cj p } , 
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which can only increase with time. The entropy S [\&] 
characterizes informational content of the market. 

From Eq. (97) we find the relation between averages: 
Vp = Qp ~ — k (see Eq. 84), which allows one to express 
parameters n and A 2 of our theory through corresponding 
moments of the trade volume V: 



2, \V\ 2 



(Aln V) 2 (Aln V) 2 



ln\V/V k \ 



(101) 



where V k ~ cr^Tk/ro is about squared bid-ask spread. 



4- Response functions 

From Langeven equation (97) we find the response of 
the mode p on volume imbalance v p (t) : 



">„(*) = / x P (t-t')v P (t')dt', 

J — OO 

Xv (t - t') = JL e -(*-*)/T,0 ( t _ t >) . (102) 




FIG. 17: The price shift AP per trade vs. transaction size 
V, for buy orders in 1996, renormalized by powers of market 
capitalization G. 91 Theoretical prediction, Eq. (105), is shown 
by solid line. Results for 1995, 1997 and 1998 are very similar. 
In Insert: normalized A T P for different r as function of V, 
V > (O) and V < (□), data from Ref. 92 . Solid lines show 
fitting by Eq. (106). 



Using Eq. (93) one can check, that the response function 
X of volatility on volume imbalance at the coarse graining 
time T r = t is determined by the sum of contributions of 
modes p > r: 



w(t)~ 
X(t-t'\r) = 



X{t-t'\r)v{t')dt', 

i 

'(t-f) 



'p=0 



X v (t - 1') 



t-t' 



(103) 



The average price shift because of a single trade at the 
time t = of the volume 

\AV\ = V k (e Av - 1) 

(— 1 is only important at small |AV| ~ Av) is determined 
by all modes with times from Tk through tq (the noise 
£ fe ~ 1 at T k ): 



AP (t) ~ o- k 



1 



(TfeAw (t) 



(104) 



CfcX(£K)r fe At>. 



The dispersion ak at time interval r = is obtained 
by averaging over fluctuations of random variables (98), 
describing variations of the liquidity of the market. In 
general, the liquidity (at physical language, susceptibil- 
ity) strongly depends on the history: small volumes can 
initiate large jumps or make almost no effect. Similar 
"aging" effect exists for spin glasses, where the suscep- 
tibility strongly depends on the history of temperature 
and magnetic field changes. 
Using Eq. (92) we get 

AP (t) = G (t) sign(AV) In (1 + | AV\ /V k ) , (105) 
G (t)~a k 6(t) r k (t + Tk)- 1 . 



In general, 1 under logarithm is out of accuracy of our 
calculations, and we hold it to reproduce expected lin- 
ear response AP ~ AV at extremely small AV. The 
result Gq (rk) — <Jk = & (i~k) extremely well supported 
by data 90 . Weak logarithmic dependence of average price 
shift AP on the trade volume AV is related to multi-time 
character of volume fluctuations, described by Eq. (92). 

The dispersion, ~ G _/3 , is inversely correlated to 
the capitalization G of the market, see Eq. (13). The 
exponent j3 ~ 0.3 (Eq. 88 at r = Tk) is smaller than the 
Gaussian value 1/2, because of hierarchical structure of 
financial markets, see Ref. 93 . It is shown in Fig. 17 that 
price impact curves for 1000 stocks are collapsed very 
well by Eq. (105) with V k - G s and S ~ 0.3 - 0.4. 

The average price shift during the time interval r > 
Tfc can be estimated considering several trades as one 
large trade of summary volume AV, and renormalizing 
minimal relaxation time — ► r in Eq. (105): 



A r P ~ asign (AV) In (1 + | AV| /V T ) 



(106) 



with V T 2 = V k 2 T/Tk- We show in Insert in Fig. 17, that 
the above dependence A r P well agrees with empirical 
data at different r. Eq. (106) can be approximated by 
power dependence AP ~ hAV 1 ^, with time, volume 
and capitalization dependent effective exponent (see In- 
sert in Fig. 17): 



dlnV 
d\nA T P 



= 1 



IAVI 



In 1 



\AV\ 



At small AV the apparent exponent v is large at small r, 
and v ~ 1 at large t (v ~ 3 for r = 5 min and v ~ 1 for 
t = 195 min, see Ref. 92 and Insert in Fig. 17). At large 
AV typical value v ~ 2, see Ref. 94 , while 1/v slowly 
decreases 95 with AV, and 96 1/v — > for very large AV. 



24 




Time (Trades) 



FIG. 18: Response function R(l,V), conditioned to a cer- 
tain volume V, as a function of dimensionless time I. Data 
for France- Telecom 90 . The thick line is the prediction of 
Eq. (107). 



In the end of this section we estimate introduced in 
Ref. 90 response function conditioned to a volume V: 



n(l,V)~Yl ^AP[(l-n)n]r,(0) 



AV=V 



\AP[(l-n) T k ]\v<JiT k )r](0)dn. 



Here |AP [(I — n) T k }\ is the value of average price shift 
at time W k because of a trade at time rvr k . It is impor- 
tant, that ln|V| at time Itu for given value of ln|y| at 
time nTk logarithmically weakly depends on time interval 
(I — n) Tk because of multi-time relaxation of this value 
(see Eq. 112 below as an example of calculation of such 
conditional average). Substituting Eqs. (105) and (96) 
we find with logarithmic accuracy 



rl 



n(i) 



±n(l)ln\V/V k \, 
1 dn ln(l + /) 



(107) 



l-n + ln"i 



l~ 



Plotted in Fig. 18 function 1Z (I) is in good agreement 
with empirical data for France- Telecom 90 with the sole 
parameter 7 = 1/5. In general, TZ(l) initially grows, 
reaching maximum at certain I* ~ e 1 / 7 , and than de- 
creasing back with dimensionless time I > I* . 

Notice, that Eqs. (105), (106) and (107) were obtained 
by pre-averaging over fluctuations of the liquidity of mar- 
ket, and can not be applied to find higher moments of 
price increments, like dispersion <t(t). Multifractality 
changes power dependences of these values: in contrast 
to the first moment, Eq. (107), when the intermittency 



effect is not important, it gives a major contribution to 
higher moments. We show in section IV C 6, that the 
liquidity fluctuations lead to strong variations of the vir- 
tual trading time, the rate of which is determined by 
local time between trades. Therefore correlation func- 
tion 7Z(l,V) with pre-averaged time intervals r k car- 
ries no information about dispersion a (r). The relation 
between 1Z (I, V) and a 2 (r) ~ Dt was used in Ref. 90 
to demonstrate a very delicate balance between liquid- 
ity takers and liquidity providers to put the market at 
the border between sub- and super-diffusive behavior. In 
section IV B 1 we show, that apparent diffusive behavior 
<j (t) ~ Dt is really a result of random trader activity 
at all time scales. 



5. Stock and news jumps 

In this section we study volatility patterns of large 
price jumps in the market. The volatility variable 
oj (t) (90) can be measured empirically as the average 
over n> 1 time intervals r of the logarithmic modulus 
of price increments: 

1 x — \ n 

u (t = - V\ In \AP k {t)\ , AP k (t) = A T P (t - kr) . 

In the limit n — > 00 uj (t) can be considered as asymp- 
totically Gaussian random variable. To prove this it is 
instructive to define generalized volatility 



v q (t) 



1 n 



(108) 



which turns to standard definition of volatility at q = 1, 
while in the limit q — > we have 



w(t) = 



dVg (t) 



dq 



(109) 



q=0 



In Appendix E we show, that at large n ^> 1 PDF 
of volatility converges to universal function, which de- 
pends only on q, exponent fi = 3 of the fat tail of PDF, 
V (AP) ~ lAPp 1 "^, and non-universal constant c > 
(will be calculated later): 



P(V q ) 



cF (c/V q) V„, 



(no) 



The maximum of this distribution is reached at the point 
V max = [c (1 + l-i/q)Y C V m . 

In Fig. 19 we show that expression (E2) of Appendix E 
for q = 1 with fi = 3 and c\ = 2/3 is in excellent agree- 
ment with known empirical data. Usually this depen- 
dence is fitted by log-normal, Eq. (E3) of Appendix E, 
or inverse gamma distributions 97 (Eq. (110) with c = 1) 
with extremely high exponent /j, = 5 — 7. Our calcula- 
tions show, that the distribution P (V q ) has the same tail 
exponent fi = 3 as PDF of price incemants. This result 
clearly demonstrates the absence of the self-averaging of 
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FIG. 19: Comparison of theoretical (with fi = 3 and c = 2/3) 
and empirical 87 volatility distributions for r = 30 min, nr = 
120 min. In insert we show the tail distribution of volatility 
of jumps and news 77 . 



price fluctuations: large variations of the coarse-grained 
"volatility" variable V q (t) (108) are induced by large 
short time jumps, the contribution of which is dominated 
even after averaging over nr — > oo time interval. 

In Insert in Fig. 19 we show the probability of large 
volatility fluctuations 77 . As one can see, the probability 
of large "stock jumps" has power tail, P (V) ~ V ~^ 
with /i — 3, while /j = 2 for "news jumps", induced 
by independent macroeconomic events. We show in sec- 
tion II C 4 that /i equals to the number of essential de- 
grees of freedom of the noise: complex uncorrelated noise 
of "news jumps" has /i = 2 components, while there is ad- 
ditional component of price at previous time interval for 
Markovian noise of stock jumps. The prediction \x = 3 is 
strongly supported by the analysis of distinct databases 
with extremely large number of records 98 for the interval 
r from a minute through several months. 

Both predictions for the tail exponent of PDF, /x = 2 
and \i = 3, are quite general. For example, they de- 
scribe two major universal classes of city grow (discov- 
ered from empirical data in Ref. 99 ), because of adding 
new street lines. The PDF describes the distribution of 
lengths of these lines. New lines are created randomly for 
cities with /j = 2, while there are strong local correlations 
in line creation for cities, characterized by the exponent 
jj, = 3. Similar Gutenberg- Richter power law describes 
earthquakes of a given strength. 

Since \xj q — > oo at q — > 0, the distribution of V q (110) in 
this limit becomes asymptotically Gaussian, and random 
variable u> (t) at large n 3> 1 has Gaussian statistics with 
the probability 

P{w}~e- ff M ) 
H {lu} = 1 Jj dtdt'io (t) to (f) G- 1 (t - f ) , 1 ' ' ' ' 

where G _1 is an inverse to the kernel G, Eq. (100) (ex- 
plicit expression for G _1 is given in section IV D). We 



checked Gaussian character of oj (t) by numerical simu- 
lations in the model of Eq. (78). Eqs. (90) and (100) 
lay in the basis of the famous Multifractal Random Walk 
model 78 ' 79 . 

Minimizing H {lu} (111) under the condition of fixed 
u> (to), we find deterministic component of uj (t) ait— to > 
t: 

lu (t) = AG (t) = lu (t ) h(t — t ) , (112) 
h(t) = eln (t / \t\), e= 1/ln (r /r) , (113) 

where we expressed the Lagrange multiplier A = 
(to) /(2A 2 ) through Lu(t ). 

Eq. (112) describes the result of trading activity, while 
news coming at time t = induce additional volatility, 
see Eq. (103): 



Lj n (t) = w r/t, 



t > T, 



(114) 



lu is the amplitude of the news jump. Combining both 
contributions, Eq. (112) and (114), we find differential 
equation for the most probable to (t) after a news jump 
du = —eLud\n(t/T) + du) n , with the solution 



Lu(t) 



1-e 



t > T. 



(115) 



The resulting volatility pattern, a (t) — aoe"^, can be 
measured by averaging over all significant news. 

Last term in Eq. (115) appears because of long mem- 
ory effect, known as "aging" effect in spin glasses. It can 
also be interpreted 77 as the reduction of measure of un- 
certainty after news, since some previously unknown in- 
formation becomes available. In Fig. 20 we show the pre- 
diction of Eq. (115) in comparison with empirical data 77 . 
Initial increase of u> (t) at t < is related to finite waiting 
time t w , see Eq. (115). 

Although stock jumps have Markovian statistics, cor- 
relations quickly decay at the coarse graining time r, 
which is slightly longer than the time of returns, and 
at larger times these events can be described by uncor- 
related random variable r/(t). The volatility, Eq. (100), 
can be rewritten through rj (t) as 100 



Alu (t) 



dt'rj (f) 



V (t)r,(t>)=\'S(t-t'). 



(116) 

Substituting 77(f) = loqS (t) in Eq. (116), we find, 
that volatility after stock jump at f = is relaxing 
as lu (t) = wor 1 ' 2 , in good agreement with empiri- 
cal observations 77 . Stock jumps can be interpreted as 
stochastic resonance of different fluctuation modes p in 
Eq. (90) because of their random concurrence. Such res- 
onance is usually happened because of delaying the jump 
until it will be anchored by jumps of larger time scales. 
The amplitude of the resulting stock jump is significantly 
increased, and slowly decays with time. The model pre- 
dicts diffusion dependence Wq ~ A 2 f tu of the amplitude 
of a jump on the waiting time t w between neighboring 
stock jumps. 
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FIG. 20: The volatility pattern before and after news jump, 
given by Eq. (115) with r = 3, to = 10 years, in comparison 
with empirical data 77 . The amplitude loo is choosed from 
the condition of best fitting. In Insert we show, that the 
probabilities of stock jump after jump, P (J, t\J, 0), decays 
to jump probability as |i| _1//2 , as predicted by our theory. 
We also show the probability to observe a jump after news, 
P (J,t\N,Q), Eq. (118), which is increased at small times and 
decreased at intermediate times. 



The central part of PDF of the volatility fluctuation 
Vi (to) — aoe w '*°\ averaged over all fluctuations, can be 
found by substituting Eq. (112) into (111): 



6. Virtual trading time 

In this section we show, that fluctuations of the liquid- 
ity of the market lead to corresponding variations of the 
virtual trading time, (t), which is proportional to the 
number of trades per given time interval. Logarithmic 
volatility to (t), Eq. (112), gives the deterministic part of 
time dependence of the amplitude at t — to > t: 

a(t)~a e"V ^ a(t Q )[(t -t Q ) /r] a , (119) 

with the "feedback parameter" 

a = -ecu (t ). (120) 

In Multifractal models 101 the (logarithmic) price is as- 
sumed to follow 

P(t)=B[G(t)}, (121) 

where B (t) is Brownian motion and (t) is the ran- 
dom trading time, which is an increasing function of t. 
Differentiating Eq. (121) over t, we can represent the in- 
crement of price in the form A T P (t) — a (i) £ (t) with 
a (t) ~ Q' (t). Substituting "classical trajectory" a(t) 
from Eq. (119), we find 

e (*) - e (t ) ~ (t - t ) 1+a , (122) 

and the mean square increment of the price is 



P(Vi. 





r a 2- 




exp 


-G(0)-_ 


~ exp 



4A 2 V «o 



(H7) 

The log-normal of this distribution is supported by nu- 
merous empirical data 97 . Comparing this expression with 
Eq. (E3) of Appendix E, we find the value of constant 
c ~ A 2 (fi + 1) /e in Eq. (110) for q = 1. 

The volatility pattern of a news jump at time t = 
followed by a stock jump at time t can be presented as 
the sum of corresponding volatility patterns, Eq. (115) 
and (112). Substituting it into Eq. (Ill) we find the 
probability of this pattern: 



P(J,t\N,0) = P n (^ )P(Vi)exp 



2A 2 



_ (H8) 

Here P n and P are probabilities of news and stock jumps, 
and function to (t) is calculated in Eq. (115). We con- 
clude, that at small time interval t there is asymmetric in- 
crease of the probability to see a jump induced by a news, 
following by the decrease of this probability at intermedi- 
ate times. Similar expression (118) with to (t) = tuot^ 1 ^ 2 
can be found for the probability P (J, t\J, 0) to find a 
jump at time t after a jump at time t — 0. These pre- 
dictions are in good agreement with empirical data, see 
Insert in Fig. 20. 



jp(t)-p(to)r;~e(t)-e(to) 

with the local Hurst exponent 

H = (l + a)/2. 



(t-tof H (123) 



(124) 



Expression (123) is valid for any to with current H(to). 

The price P (t) experiences different types of fractional 
Brownian motion in time intervals Ati with different 
feedback index a, which randomly change each other, 
see Fig. 21. The case H ~ 1/2 (a ~ 0) describes usual 
Brownian motion. A Hurst exponent value < H < 1/2 
(a < 0) will exist for a time series with sub-diffusive 
(anti-persistent) behavior. A Hurst exponent value from 
the interval 1/2 < H < 1 (a > 0) indicates super- 
diffusive (persistent) behavior. H ^ 1/2 can be inter- 
preted as the result of local unbalance between the com- 
peting liquidity providers and liquidity takers 90 . 

The fractal dimension of the fractional Brownian mo- 
tion is D ext — 2 — H . At large H the motion becomes 
more regular (D ext — > 1), with large up- and downturns, 
while at small H it quickly fluctuates, trying to covers 
the whole plane (D ex t — > 2). Therefore, establishing of a 
super- and sub-diffusive behavior leads to significant sup- 
pression/creation of short-time fluctuations, see Fig. 21. 
This effect was really observed 102 , and may be used as an 
indicator of establishment of large volatility at long time 
scales, which is hard to detect at short time intervals. 
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FIG. 21: Price fluctuations and corresponding virtual trading 
time Q (t) as functions of real time t for Brownian motion 
(a — 0), sub- (a < 0) and super-diffusive behavior (a > 0). 



FIG. 22: Daily series of Exxon Mobil Corporation and corre- 
sponding \i (t) 102 . Probability distribution of ai — hiao with 
hi = 0.9 at Ati = 16 days is shown in Insert. 



7. Brownian motion, sub- and super- diffusion 

Switching of fluctuation regimes between Brownian 
motion (a ~ 0), sub- (a < 0) and super-diffusion (a > 0) 
is happened randomly at "frustration times" (with equal 
probability of different choices) with the probability (see 
Eq. (117)) 



p(a) 



I 



exp 



2o\ 



al = 2eX 2 . (125) 



Because of multifractality, such change of fluctuation 
regimes is happened at all time scales r. 
We also define a multivariate PDF 



transition to a new state is happened in average when 
two these amplitudes become of the same order: 



Ah 



rj- 

V T 



i/ym 



4e 2 A 



2 ■ 



(130) 



It is surprising, that the length of the time interval At\ 
grows with the rise of |ct | (although the probability of 
large initial |arj|, Eq. (125), is small). Such counter- 
intuitive behavior is related to the absence of any "restor- 
ing force" to a = 0. 

The characteristic time of sub- and super-diffusion 
behavior can be roughly estimated substituting in the 
above expression the most probable value a\ ~ a\ from 



p {a , ■ ■ ■ , a fc ) = ]J ;=0 6 \ a l + euj ( 126 ) E q- ( 125 )> g ivin g * - \ m ( r o/r). For r = 1 day, r ~ 10 



which determines information entropy conditional to the 
set of indexes a>o, ■ ■ ■ ,ak- 



S (a , ■■■ ,a k ) ~ Inp(a , • • • ,a k ). 



(127) 



The entropy, Eqs. (127) and (125), is maximal for Brow- 
nian motion, a = 0, sub- and super-diffusion lower the 
entropy production, since the market behavior is more 
predictable in these regimes. 

Conditional dynamics of mode switching can be de- 
scribed by the probability p(ai|ao) = p{&o,(Xi) /p(ao) 
to find given value of the index a.\ at time t\ = to + Aii 
under the condition that it was arj at previous time to- 
Calculating the average (126) at k = 1 with Gaussian 
distribution function, Eq. (Ill), we find: 



p{a 1 \a ) 



1 



2ir<r 



■ exp 



(qi - Oi) 
2aj 



(128) 



ai = (Ati) , tr? = trg [l - /i 2 (Air)] (129) 

Conditional average dt\ decreases with time A^, while 
the conditional dispersion o\ grows with this time. The 



days and A 2 ~ 0.1 this gives At\ about a month, in agree- 
ment with empirical observations, see Fig 22 and Ref. 102 . 
The effect of sub- and super-diffusion switching can be 
hardly detectable at small time intervals r, because of 
small a (see Eq. (120)), but it is well pronounced for 
daily time intervals r. 

The feedback index a is related to "variation index" 
/i = (1 — a) /2, introduced in Ref. 102 from the fractal 
analysis of empirical data (local extension of the R/S 
analysis 103 of the Hurst exponent H ) . In Insert in Fig. 22 
we plot estimated probability distribution of ct\ — h\ao 
for Aii = 16 days, which is proportional to the condi- 
tional probability (128). The value hi is chosen to get 
a maximum at the beginning, and it decreases with the 
rise of the time interval Ati from 0.99 at At\ = 1 to 
0.9 at Ati = 16 days. The dispersion of this distribu- 
tion (To — 0.2, and from Eq. (125) we get reasonable 
estimation A 2 ~ 0.08. From second Eq. (129) we find 
the dispersion o\ = a$ yl — h\ ~ 0.07, close to observed 
value. 

We can also calculate the probability to find the feed- 
back index a k at time t k under the condition that it was 
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a/c-i at time tk-i, Uk-2 at time tk-2, and so on: 

,a k -i) (131) 



p(a fe |a , • • • ,afe-i) 
= p(ao,--- ,«fe)/p(ao, 



27TCr fe 



exp 



(a fc - afe) 
2^ 



The logarithm of this probability determines the entropy 
lowering because of the knowledge about previous events 
ao, ■ ■ ■ ,a k . The conditional average a k corresponds to 
the maximum of entropy production, Eq. (127). It has 
the meaning of average response of a on previous values 
a , - • • ,a k -\: 



fc-i 



<3fe 



i=0 



-dki/akk, (132) 



afej are adjoints of the matrix h with elements ft. (^ — tj). 
Conditional dispersion 



a\ = <rl dct h/a kk 



(133) 



estimates the accuracy of the prediction (132). Probabili- 
ties (131) depend not only on index a k -\ at previous time 
t k -i, but also on all ao, • • • , Qfe-i _ the random process 
is not Markovian. As the consequence, the probability to 
have the same value of all three indexes <xi = a\ = a 
(continuation of a sub- and super-diffusive regimes) grows 
with increase of the initial time interval Aii . 



8. Fluctuation corrections 

The "classical trajectory" (119) predicts that the feed- 
back index a (120) can be changed only because of fluctu- 
ations. In this section we demonstrate, that fluctuations 
lead to an additive shift of a in the super-diffusion di- 
rection. Calculating the average of V q (t) with Gaussian 
weight (111) under the condition lo (0) = cjq, we find: 



P(t) = \ 2 [l + h(t)}/2, 



qa+q 2 0(t) 



(134) 



where h (t) is defined in Eq. (113). The /3-term describes 
deviation from fractional Brownian motion because of 
multifractal behavior. At q = 2 this expression can be in- 
terpreted as the response on "endogenous event" lo (0) 100 , 
while at q = 1 it gives fluctuation correction to the feed- 
back index a. 

Empirical study of year correlations shows 104 , that 
financial market is really "locked" in sub- and super- 
diffusive or Brownian motion states at extremely long 
periods (conventions can persist up to 30 years). The 
change in convention can be rather smooth, like during 
the second part of the century, or occur suddenly, trig- 
gered by an extreme event, like it did after 1929. From 
the data, it was also observed a systematic bias towards 
the persistent following convention. 



D. Universality of fluctuations 

Non-universal properties of market at trading times < 
r k can only be described by models of agent-based strate- 
gies. In this section we consider only universal properties 
of price fluctuations at time scales r > r k ~ 1 min, see 
Ref. 105 . At qualitative level the presence of universality 
is known for a long time as "stylized facts" 40 . We show, 
that our approach captures such stylized facts and gives 
explanation for many others empirical observations. 

The universality is related to self-similiarity of price 
fluctuations at different time scales 79 ' 106 ' 107 : the change 
of time interval r corresponds to the change of char- 
acteristic scale along the hierarchical tree in Fig. 14. 
We demonstrate, that resulting time series have com- 
plex non-periodic behavior with chaotic changes of usual 
Brownian motion, sub- and super-diffusion, reflecting 
cyclic dynamics of the market. We show, that the im- 
pact function of the market logarithmically depends on 
volume imbalance. 

Fluctuations on financial market have unexpected 
physical interpretation, reflecting the unified nature of 
physics. The effective Hamiltonian (111) can be rewrit- 
ten as 



H{w} 



JL 

2ir 



u(t)-w (*') 

7~t< 



dtdt'. 



This expression describes diffusion of quantum Brownian 
particle with the coordinate u> (t) and the coefficient of 
linear friction r] = 1/ (27rA 2 ). 

A microscopic model of quantum diffusion is based 
on coupling to a termostat - the reservoir of harmonic 
oscillators 108 , presenting the "army" of traders in the 
case of the market. The resulting dynamics is intrinsi- 
cally non-Markovian in that the evolution depends on his- 
tory rather than just on present state 109 . Brownian parti- 
cle can respond to a very wide range of reservoir frequen- 
cies, and this is the origin of time-irreversive behavior and 
slow relaxation after fluctuation cast, see Eq. (112). The 
production of information entropy (see Eq. (127) and Ap- 
pendix A) is related to enviroment-induced decoherence 
of the quantum particle 110 , and it is at the peak of many 
recent studies. 



V. CONCLUSION 

Many concepts of equilibrium macroeconomic (re- 
sources, unemployers, different firm dynamics at small- 
and long- time horizons, taxes and so on) naturally 
enter into proposed coalescent theory, which integrates 
both physical and economic concepts of essentially non- 
equilibrium market in one unique approach. We also de- 
veloped new approach to study fluctuations on the mar- 
ket, well describing empirical data of both firm grow rates 
and price increments on financial markets. We propose 
the set of Langeven equations, describing multi-time dy- 
namics of price and volume fluctuations at different time 
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scales on the market. Using these equations, we de- 
rived analytically equations of multifractal random walk 
model. 

In the end we discuss physical meaning of our theory: 

a) What is physics of extreme events on the market 
(problem of fat tails)? 

There are two sources of fluctuations: macroeconomic 
events - news, and traders activity because of uncertainty 
of equilibrium prices at different time scales, which gen- 
erate two types of large price jumps: News jumps are 
created by the inflow of news, while stock jumps are gen- 
erated during random concurrence of different fluctuation 
modes. We show, that the decay of volatility observed af- 
ter news jumps is related to the effect, similar to "aging" 
effect in spin glasses. 

Fluctuations on the market are characterized by the 
normalized noise £ and its amplitude a (volatility) . New 
key idea of our approach is that £ and a are independent 
complex random variables, separated on time scale: the 
noise is generated by hot degrees of freedom on times 
small with respect to observation time interval r, while 
evolution of the amplitude is determined by cold degrees 
of freedom on times large with respect to r. 

In stochastic volatility and multifractal models jumps 
are predicted as the result of volatility fluctuations, and 
characterized by large non-universal tail exponent /x >• 1, 
while the noise is assumed to be Gaussian uncorrelated 
random variable. In fact, the noise is strongly correlated, 
and can experience large non Gaussian jumps. We cal- 
culate the contribution of such jumps to PDF and show, 
that the distribution of stock jumps is characterized by 
the tail exponent fj, = 3, while the distribution of news 
jumps has tail exponent fi = 2. The exponent fi remains 
stable with the rise of r (recall, that Levy distribution 
with fj, > 2 is unstable). 

b) Why market dynamics is so complex: can it be de- 
scribed by simple Markovian or Gaussian processes? 

We show, that local equilibriums on the market are 
self-organized in the hierarchical tree, according to their 
relaxation times. The amplitude a of the noise at given 
time scale is determined by cumulative signal from all 
"parent" time scales, and its dynamics is complex multi- 
fractal process. But the information about the amplitude 
can be "erased" from time series considering only signs of 
price fluctuations. The resulting Markovian process de- 
scribes propagation of positive and negative signals, and 



determines conditional double dynamics of the market. 

Typical noise £ and amplitude a are determined by sig- 
nals from large number of, respectively, short and long 
(with respect to r) time scales, and they have asymptot- 
ically Gaussian statistics. We propose and solve Dou- 
ble Gaussian model of market fluctuations, and show 
good agreement with empirical data for different groups 
of stocks. 

c) What physics stands behind "random trading time" 
in Multifractal models 101 ? 

At large time intervals the price randomly cycles 
between Brownian motion, sub- and super-diffusive 
regimes, which change each other because of liquidity 
fluctuations. The virtual trading time is proportional to 
the real time for Brownian motion and experiences time 
shifts in sub- and super-diffusive regimes. The theory 
predicts systematic bias to persistent behavior, observed 
for many markets and exchange rates. 

d) And finally, can price behavior be described by uni- 
versal physical lows or it is dictated only by the zoo of 
microstructures of markets (see Refs. 111,112 )? 

The universality of price fluctuations on financial mar- 
kets was demonstrated at time scales from a minute 
to tenths years in many studies, see for example, 
Refs. 28 - 36 - 37 . We show, that it is related to the self- 
similarity of the underlying hierarchical tree of ampli- 
tudes, see Fig. 14 (we do not give here lists of all stocks, 
used for comparison with our theory, since they are shown 
in corresponding references). In contrast, statistics of 
trades and volumes is not universal, and strongly de- 
pends on details of market microstructure. 

Our theory can also be used to study other time series, 
such as variations of cloudiness, temperature, earthquake 
frequencies, rate of traffic flow and so on. It looks attrac- 
tive to apply analytical approach of this paper for the 
description of social processes, which are driven by frus- 
trations at turning points of the mankind history. Events 
between these points support the social activity, but do 
not change the state of the society. 
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APPENDIX A: ENTROPY FORMULATION 

In order to reveal the economic meaning of Master 
equation (10), it is convenient to rewrite it in the form 



dC 



(Al) 



where p c (G) and pd (G) are probabilities of job creation 
and destruction per unit time in the firm of G people. In 
the absence of any external supply, U = 0, the probabil- 
ity of job creation is zero. In the main order in "concen- 
tration" U (we use physical term to emphasize the anal- 
ogy with coalescence) p c (G) is proportional to U, while 
the probability of job destruction p d (G) is determined 
mainly by internal firm structure, and do not depend of 
U. Comparing Eqs. (Al) and (10) we find explicit ex- 
pressions for these probabilities for our model 



Pc (G) = qUG, Pd (G) =qU*G+ P G 



1-/3 



(A2) 



We define the "entropy" S (G) of the firm of size G 
as logarithm of equilibrium distribution function of firms 
over their sizes, f eq (G). Since equilibrium values do not 
depend of a way how the system is assembling, consider 
the process when the size G is varying by one. In this case 
feq (G) is determined by the detailed balance condition, 
Pc (G) f eq (G) = Pd (G + 1) f eq (G + 1). The solution of 
this equation atG>l relates firm entropy with proba- 
bilities of job creation and destruction: 



f eq (G) = e s ( G \ S(G) 



dG In 



Pc(G) 
Pd (G) ■ 



(A3) 



As one can naively expect, the entropy of the firm in- 
creases with the rise of the probability to get a job and 
decreases with the rise of the probability to loose it. In 
the case of overheated market U < U* the entropy S (G) 
monotonically grows with the firm size G, while in the 
"supersaturated" case U > L 7 * it initially decreases with 
G, reaching its minimum for firms of critical size, G = G c , 
Eq. (14). 

Calculating this integral (A3) with functions (A2), we 
find 



S(G) ~ uG - Gin \u./Uo + (G/e) 



+ const, 



where Uo = p/g and u = \n(U/Uo) has the meaning of 
chemical potential. The entropy of the whole market 

S = -Uln-^-+ f dGS(G)f(G,t)-u(Q-U) (A4) 
eUo J 

is the sum of the entropy of the "ideal gas" of unemploy- 
ments, the entropy of all firms and the term <~ u, which 
takes into account the supply of external resources (1). 
By analogy with thermodynamics, there is maximum 
principle for the entropy: maximizing it with respect to U 
we reproduce the balance condition (1). Using Eqs. (A3) 
- (A4) one can check that in the case when Q (t) is not 
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(quickly) decreasing function, the market entropy always 
increases with time 



dS/dt > 0. 



(A5) 



Since the variation of entropy AS is opposite to the vari- 
ation of information, AI = —AS, Eq. (A5) means, that 
the activity of the market leads to "erasing" of initial in- 
formation - the effect, well known for some "laundering" 
schemes. 



APPENDIX B: SOLUTION OF COALESCENCE 
EQUATIONS 

To find PDF of coalescent model (10) we use the 
method of Ref. 113 . We define dimensionless time r and 
introduce the function u (t): 

t = In [G c (t) /G c (to)] , u(t) = G (t) /G c (t) , 

where to is coalescent time. In new variables the Master 
equation (10) takes the form 



du/dr = v (u) = 7 (t) (u - u 1 -?) - u, (Bl) 



where 



dt 



7(T)= W 



(B2) 



The balance equation (5) can only be satisfied if the plot 
of function v (u) lays below the axis u, and touches it at 
one point u = u$. Such locking point u = uo, 7 = 7 , for 
Eq. (Bl) is determined by equations 

v (uq) — 0, dv (uo) /duo = 0, d 2 v (uq) /duQ < 0, 

and we find that uq — > 00 and j = 1. 

From Eq. (B2) we get that the critical size G c (t) grows 
and "supersaturation" A (t) (we use physical terms here) 
decreases with time as 



1//3 



A(t) = J t . (B3) 



PDF of firms can be rewritten through PDF of vari- 
ables u and t: / (G, t) = <p (u, r) /G c (t), and neglecting 
the diffusion inflow of new firms we find 



dip d , , . 



where the velocity v (u) = du/dr 
Eq. (Bl) with 7 = 7o . 

General solution of Eq. (B4) is 

<P{u,t) = t/ _1 X [t - t {u)\ , 



(B4) 

1-/3 is given by 



t(u) = -u> 3 /(3 (B5) 

with arbitrary function x ( T ) • To find \ ( T ) substitute 
this expression into the balance equation (1) and (5) with 
Q(t) » U(t): 



(*o)e 



r(/3m-l) _ 



duuip (u, t) . 



This condition can be satisfied only if the function % 
has the form X [t - t (u)] = Ae^" 1 - 1 ^-^. Substitut- 
ing this function into Eq. (B5), we find 



with 



(«, t) = A (1 - (3m)- 1 e^ m - 1)T dF (u) /du, (B7) 



(B8) 



APPENDIX C: MACROECONOMIC 
INTERPRETATION 

To get better understanding of coalescent model, con- 
sider its microeconomic interpretation. Optimal firm 
size, G = G c , is determined from the maximum of the 
profit function 



tt (G) - Py (G) - wG, 



(CI) 



where y (G) is the number of units produced by firm 
of G peoples, P is the price of one unit, and w is the 
average wage per one man. The technology is usu- 
ally characterized by the standard Cobb-Douglas func- 
tion y (G) = K ri G 1 ^ n , where K is the firm capital and 
the exponent 77 > 0. Both the price P and the capi- 
tal K are reduced to initial time. Maximizing the profit 
function (CI) we get G c ~ w" 1 ^. 

Variation of wages w (t) with the time is determined 
by the Fillips low 114 : 



1 dw (t) 
w (t) dt 



~a[U*-U (t)] = -oA (f) . 



(C2) 



(B6) 



with positive constant a > 0. Substituting expres- 
sion (B3) for A (t) in Eq. (C2), we find its solution 
w (t) = const x t _< » with £ = a/ (J3q). Substituting this 
function into G c <~ we get the optimal firm size 

G c <~ i^/ r '. Comparing this dependence with Eqs. (14) 
and (B3), we find the Fillips parameter in Eq. (C2): 
a = nq. We conclude, that coalescent approach is consis- 
tent with the maximum profit principle and the Fillips 
low. 

Notice that while the parameter q = a/rj of Master 
equation (Al), (A2) is determined by technology (77) and 
market structure (o), economic analysis do not impose 
any restrictions on the second parameter p of Master 
equation. Therefore, p can depend on management abil- 
ity of firm head, relations between firm staff, industry 
shocks and so on, and can experience strong random 
fluctuations Ap. This observation explains the empirical 
fact, that there is much more variance in job destruc- 
tion than in job creation time series 115 (as was noted by 
Lev Tolstoy: all fortunate families are happy alike - each 
unfortunate family is unhappy in own way). 
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APPENDIX D: PDF OF DOUBLE GAUSSIAN 
MODEL 

It is convenient in Eq. (62) to use instead of {a°} Gaus- 
sian random variables {on}: 



a, = CiOLi 



such as af = 1, (ai, c* 2 ) = v. Averaging over fluctuating 
Gaussian variables £° and on, we get general expression 
for inverse Fourier component of PDF: 

G~ l (k,p) = 1 + k 2 a u /2 + p 2 a 22 /2 + kpa 12 + 



with 



(l - v 2 ) {xiik 2 + K 2 2P 2 + Knkpy 



Eq. (D2) can be simplified if we note, that correlations of 
price increments, Eq. (Dl), are always very small, \e\ -C 
1. In polar coordinates K = |K|cosy>, P = |K|sin<^ we 
have in the main order in e 



\KP-t{K,P)] 
1 

4 



K 



- 4 K| 4 sin 2 [2(^- 



sin (2tp) 

COS (2(f)))] = 



£COS (2(f))Y 

(K'P'f , 



where we introduced new orthogonal rotated coor- 
dinate system K ~ K' - (e/2) P' cos (20) , P ~ 
P' + (e/2) if' cos (20). Changing integration variables 
— > (K',P') in the integral (48) in the main or- 
der in the small parameter e we get Eq. (63), where an- 
gles + and 0_ are defined by 0_ = — e cos 2 0, + = 
4> + e sin 2 0, and functions "P; (a;) are determined by 
Eqs. (58) and (67). 



(Til =Cj_ + C 12 + C n + 2l/CnCl2, 

^22 = C 2 + C 21 + C 22 + 2I/C22C21, 

(T12 = C1C21 + C 2 Cl2 + f (ciC 22 + C 2 Cn) , 

^11 = C1C12, >* 22 = C2C2I, 

X12 = ClC 2 - CnC 22 + C12C21. 

An important relation between elements of matrixes a 
and x follows from the condition of stationarity of PDF, 
which leads to physical constraint G (k, 0) = G (0, k) on 
Fourier components of univariate PDFs of A T P (t) and 
A T P (t + r), and we get crn = (T 22 = cr 2 ,^^ 
Here a is the dispersion of price fluctuations, 



^22- 



(A r P 2 (i)) = (A r P 2 (i + r)) = a 2 



(A T P (t) A T P (t + t)) = a 12 = Ed 



(Dl) 



We conclude, that in addition to a and e, our model is 
characterized by dimcnsionlcss constant ^ and the angle 



G- 1 [k,p) = 1 + [{k 2 + p 2 ) 12 + ekp] a 2 + (l - i/ 2 ) x 
' — — sin (20) /2 - fcp cos (20) . 



T 



(D2) 

The quadratic part of this expression can be diagonal- 
ized by changing variables, K = kcosip ~ ps'mip' , P = 
fcsin'0 + pcosip', with ip' = cfr — ^ arcsine, -0 = + 
i arcsine. In new variables Eq. (D2) takes the form 



G- 1 (k, P ) = l 
1 



P 2 + 



j/ 2 a 4 



2\2 



[KP-r(K,P)} 2 , 



with 



r (if, P) = £ cos (20) (K 2 + P 2 ) + l - (1 - ^/l _ ^) 
x sin (20) [2APsin(20) + (A' 2 - P 2 ) cos (20)] . 



APPENDIX E: PDF OF VOLATILITY 
FLUCTUATIONS 

PDF of the volatility variable Vi (t) (108) at q = 1 
can be expressed through the n-point PDF P{APfe} of 
variables Apj, 



Pn {Vl) 



f[ jdAP k S (v, - I l AP *l) V { Ap fe> 

(El) 

Asymptotes of P„ (Vi) can be found both for Vi -C a 
and for Vi > a. At Vi < only small |AP fe | < er 
contribute to the integral (El) and we have P„ (Vi) <~ 
Vi" -1 . In the opposite case of large Vi » a the integral 
is dominated by the power tail of PDF P{Apj} with 
typical |APfe| ~ Vi ^> a. From dimension consideration 
we find for such AP& that P{APfc} ~ Vi~"~ M , where 
\i is the exponent of one-point PDF, Eq. (36), and the 
integral (El) is estimated as P n (Vi) ~ Vi~ 
Both these limits are matched by the function 



Pn(V 1 ) = V - 1 Nv 1 f(V 1 /V ), 

f(z) = z~ 1 [z- n/s + z» /s y 



(E2) 



with Vb ~ u. The dependence of a new parameter s > 
on n will be found later from the condition that at 
large n the distribution P n (Vi) should not depend of n. 
Momentums of this distribution (Vf) = V^Nk/No are 
determined by normalization integrals, 

00 

Nk = J z k f (z) dz = mB [m (n + k) ,m (p - k)] , 


m = s/ (n + /j.) , 

where B is the Beta-function. 

The function / (z) (E2) reaches its maximum at z max = 
[n — l) m / (p + l) m . The central part of the distribution 
is obtained by expanding the probability P„ (Vi) over 
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ln(Vi/V max ) near its maximum at V max = z maK V , and 
it has log-normal form: 

lnP n (V 1 ) = const-l— - (f, + 1) (in -p-) . (E3) 
Since the distribution P n (Vi) (E3) should not depend on 



n at large n, we find s = c (n — 1) with certain constant 
c. In the limit n — > oo P n (Vi) becomes universal func- 
tion of Vi/Vmax, Eq. (110) with q = 1. Repeating our 
calculations for general g > 0, we find that it is given 
by the substitution fi — > /i/g and n — > n/g in the above 
expressions. 



